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Abstract. Partial localization is the phenomenon of self-aggregation of mass 
into high-density structures that are thin in one direction and extended in the 
others. We give a detailed study of an energy functional that arises in a simpli- 
fied model for lipid bilayer membranes. We demonstrate that this functional, 
defined on a class of two-dimensional spatial mass densities, exhibits partial 
localization and displays the 'solid-like' behavior of cell membranes. 

Specifically, we show that density fields of moderate energy are partially 
localized, i.e. resemble thin structures. Deviation from a specific uniform 
thickness, creation of 'ends', and the bending of such structures all carry an 
energy penalty, of different orders in terms of the thickness of the structure. 

These findings are made precise in a Gamma-convergence result. We prove 
that a rescaled version of the energy functional converges in the zero-thickness 
limit to a functional that is defined on a class of planar curves. Finiteness of 
the limit enforces both optimal thickness and non-fracture; if these conditions 
are met, then the limit value is given by the classical Elastica (bending) energy 
of the curve. 



1. Introduction 

In this paper we study the asymptotic expansion as e ^ of the functional 

j e [ \\/u\ + -di{u, v) if {u, v) e /Ce, 

Te{u,v) < ^ e (1-1) 

I oo otherwise. 

Here •) is the Monge-Kantorovich distance (Definition 3.1), and 
/Ce := |(u,w) e BV{R^;{0,e-^}) x (R^; {0, e"^}) : 



J u = J V = M, uv = a.e.j . (1.2) 



How this singular-perturbation problem is related to the terms in the title — partial 
localization, lipid bilayers, and the Elastica functional — we explain in the rest of 
this introductory section. 

1.1. Lipid Bilayers. Lipid bilayers, biological membranes, are the living cell's 
main separating structure. They shield the interior of the cell from the outside, 
and their mechanical properties determine a large part of the interior organization 
of living cells. The main component is a lipid molecule (Fig. 1) which consists 
of a head and two tails. The head is usually charged, and therefore hydrophilic, 
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Figure 1 . Lipid molecules aggregate into macroscopically surface- 
like structures 

while the tails are hydrophobic. This difference in water affinity causes lipids to 
aggregate, and in the biological setting the lipids are typically found in a bilayer 
structure as shown in Fig. 1. In such a structure the energetically unfavorable 
tail-water interactions are avoided by grouping the tails together in a water-free 
zone, shielded from the water by the heads. Note that, despite the structured 
appearance of the bilayer, there is no covalcnt bonding between any two lipids; the 
bilayer structure is entirely the result of the hydrophobic effect. 

As planar structures, lipid bilayers have a remarkable combination of solid-like 
and liquid-like properties. On one hand they resist various types of deformation, 
such as extension, bending, and fracture, much in the way a sheet of rubber or 
another elastic material does: for instance, in order to stretch the bilayer, a tensile 
force has to be applied, and when this force is removed the structure will again 
relax to its original length. With respect to in-plane rearrangements of the lipid 
molecules, however, the material behaviour is viscous, and there is no penalty to 
large in-plane deformations. 

While this interesting combination of properties is clearly related to the chemical 
makeup of the lipids — most importantly, the hydrophobic character of the tails — 
quantitative and detailed understanding of the phenomenon is still lacking. Here 
we focus on a simple question that has already been alluded to above: how can 
we understand the stability of these planar structures, and their pseudo-solid be- 
haviour, if they are constructed from independent, non-bound molecules? This is 
the main question behind the analysis of this paper. 

1.2. Partial localization. The self- aggregation of lipids into bilayers is funda- 
mentally different from the self-aggregation of 'simple' hydrophobic compounds in 
water. Lipid bilayers are thin structures, in the sense that there is a separation of 
length scales: the thickness of a lipid bilayer is fixed to approximately two lipid 
lengths, while the in-plane spatial extent is only limited by the surroundings. By 
contrast, 'simple' hydrophobic compounds aggregate in water to form drops that 
lack the small intrinsic length scale of lipid bikwTs. 

We use the term partial localization for the self-aggregation into structures that 
are thin in one or more directions and 'large' in others. The word localization is 
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taken from the literature of reaction-diffusion equations, in which locahzed solutions 
are those that are concentrated, in a well-defined way, in a small neighbourhood of 
a point. By partial localization we refer to localization to the neighbourhood of a 
set that has intermediate dimension, i.e. dimension larger than zero (a point) but 
smaller than that of the ambient space. 

1.3. Energy on a mesoscale: The functional J^^. The functional J^^ that 

is defined above is the result of an attempt to capture enough of the essence of 
lipid bilayers to address the issues of stability and pseudo-solid behaviour while 
keeping the description as simple as possible. The derivation of !Fe, which is given 
in Appendix A, starts with a simple two-bead model of lipid molecules, inspired 
by well-known models of block copolymers [35, 10, 29], and proceeds to reduce 
complexity through a number of sometimes radical simplifications. Despite these 
simplifications the physical origin of the various elements of jFg remains identifiable: 

• The functions u and v represent densities of the (hydrophobic) tail and 
(hydrophilic) headb eads; 

• The term J |Vu|, coupled with the restriction to functions u and v that 
take only two values, and have disjoint support, represents an interfacial 
energy that arises from the hydrophobic effect; 

• The Mongc-Kantorovich distance di{u,v) between u and w is a weak rem- 
nant of the covalcnt bonding between head and tail particles. 

The parameter e appears in (1.2) as a density scaling. We prove that e in fact is 
related to the thickness of structures with moderate energy. 

We show in this paper that the functional J7> favours partial localization, i.e. 
that pairs (w, v) for which J-"e(w, v) is not too large are necessarily partially localized. 
We will also show that displays three additional properties that were already 
mentioned: resistance to stretching, to fracture, and to bending. To describe the 
last of these in more detail we briefiy return to biophysics. 

1.4. Energy on the macroscale: The Helfrich Hamiltonian and the EIeis- 
tica functional. Intrigued by the shape of red blood cells Canham, Helfrich, 
and Evans pioneered the modelling of lipid bilayer vesicles by energy meth- 
ods [17, 31, 24]. The name of Helfrich is now associated with a surface energy 
for closed vesicles, represented by a smooth boundaryless surface S, of the form 



Here k, k, and Hq are constants, H and K are the (scalar) total and Gaussian 
curvature, and is the two-dimensional Hausdorff measure. This energy func- 
tional, and many generalizations in the same vein, have been remarkably successful 
in describing the wide variety of vesicle shapes that are observed in laboratory 
experiments [45]. If we add an assumption of symmetry, in which both 'sides' of 
the surface have the same properties, the 'spontaneous curvature' Hq vanishes. For 
closed surfaces in the same homotopy class, by the Gaufi-Bonnet Theorem, -BHeifrich 
is essentially equal to the Willmore functional [51] 



which arises in a wide variety of situations other than that mentioned here. 




(1.3) 




(1.4) 
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In this paper we discuss the above mentioned bilaycr models in two space dimen- 
sions, i.e. we consider functions whose support resembles fattened curves in M^; in 
the limit e ^ the thickness of these curves tends to zero. If we consider the limit 
curves to be two-dimensional restrictions of cylindrical surfaces, we observe that 
the Gaussian curvature K vanishes. If we also assume that the bilayer is symmetric 
then the Helfrich energy reduces to the Elastica functional W, the classical bending 
energy of the curve, 

Wil) = I J^^^' dn\ (1.5) 

where k is the curvature of the curve. This functional has a long history going 
back at least to Jakob Bernoulli; critical points of this energy are known as Euler 
elastica, see [49] for a historical review. 

1.5. The main result: the singular limit e — > 0. In the main result of this 

paper we study the limit e ^ 0, and connect all the aspects that were discussed 
above. The result revolves around a rescaled functional 

where M is the mass of u and v (see (1.2)). 

Take any sequence (ueji'e)e>o C /Cg. If the rescaled energy Gs{ue,Vs) is bounded, 
then 

(1) the sequence converges as measures; 

(2) the limit measure can be represented as a collection of curves (which may 
overlap); 

(3) the support of is of 'thickness' approximately 2e; 

(4) in particular the total length of the curves is M/2; 

(5) each curve in the collection is closed. 

Boundedness of Ge along the sequence therefore implies partial localization: the sup- 
port of Ue resembles a tubular e-neighbourhood of a curve. In addition, stretching 
and fracture are also represented in this result: stretching corresponds to deviation 
in total length from the optimal value M/2, and fracture to creation of non-closed 
curves. Neither is possible for a sequence of bounded 0^; both are penalized in Ge 
at an order larger than 0(1), or equivalently, at an order larger than 0{e^) in J^^. 
The resistance to bending arises as the Gamma-limit of the functional Gs itself: 

(1) For any sequence {ue,Ve)e>o ClCe, 

W{T) := V W(7) < limmiGe{ue,Ve), 

' £— »0 

where F is the collection of limit curves introduced above, and W{'y) is the 
curve bending energy defined in (1.5); 

(2) For any given collection of curves T, there exists a sequence {ug, U£)e>o C /Ce 
such that >V(r) = Umg^o Ge{Ue, Ve). 

The details of this result are given as Theorem 4.1. 

Summarizing, only partially localized structures can have moderate energy, and 
such thin structures also display resistance to bending, stretching, and fracture. 

So far we have described the result that we prove in this paper. The proof, 
however, also suggests a more refined result, that we present here as a conjecture: 
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if a sequence {ue,Ve)s>o C ICg converges to a collection of curves F, then 

:Fe{Ue, Ve) > 2M + &M{r) + e^{T) + ^(^) + ^^^^^^ ^^'^^ 

In this development, 

• the function £M(r) is a measure of the deviation of ^^^j-^cngth{-f) from 
length M/2: £M(r) > 0, and £M(r) = only if this sum equals M/2; 

• 5B(r) is zero if all curves 7 G F arc closed, and strictly positive (a positive 
constant times the number of cndpoints) if at least one 7 € F is open. 

As for the Gamma-convergence above, the inequality in this statement should be 
interpreted as a 'best-possible' behaviour; while for given limit curves 7 sequences 
exist for which this inequality is an equality, it is also possible to construct sequences 
along which the difference between the left- and right-hand sides is unbounded. 



1.6. Strategy of the analysis and overview of the paper. At first glance it 
is not obvious that the functionals encode bending stiffness effects. Curvature 
terms in fact appear at order of an asymptotic development, which asks for a 
careful analysis. Most information is hidden in the nonlocal Monge-Kantorovich 
distance term. The special property that makes our analysis work is that this 
distance decouples into one-dimensional problems. This reduction is the major 
step in our analysis and the connection to the optimal mass transport problem is the 
most important tool here. Whereas the idea becomes very clear in simple situations 
which resemble the lim-sup construction in the Gamma-convergence proof, some 
effort is required to prove the lim-inf estimate in the general case, where certain 
'defects' such as high curvature and vanishing thickness of approximating structures 
have to be controlled simultaneously. 

In the next section we motivate our choice of scaling in the functional and 
illustrate the behaviour of and Ge in two particular examples. In Section 3 we 
give definitions of the Monge-Kantorovich distance and 'systems of curves' which 
describe the class of limit structures. A precise statement of our results is given 
in Section 4. We briefly review the optimal mass transport problem in Section 5. 
In Section 6 we provide a short guide to the structure of the proof of our main 
result. Theorem 4.1, and the proof itself is given in Sections 7 and 8. We conclude 
in Section 9 with some discussions of generalisations and the wider implications of 
this result. 
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Summary of notation 





functionals describing the mesoscale energy 


(1.1) 




rescaled functionals 


(4.1) 




domain of J^e , 


(1.2) 


W(7) 


bending energy of a curve 7 


(1.5) 


w(r) 


generalized bending energy of a system of curves T 


(3.9) 




Monge-Kantorovich distance 


(3.2) 


dpi-.-) 


p-Wassorstoin distance 


(3.5) 


supp(r) 


support of a system of curves F 


Def. 3.3 




multiplicity of a system of curves T 


Def. 3.3 


|r| 


total mass of a system of curves F 


Def 3.3 


LiPi(»') 


Lipschitz continuous functions with Lipschitz constant 1 


(5.3) 


T,£ 


transport set and set of endpoints of rays 


Def 5.4 


E 


{s : 7(5) lies inside a transport ray} 


Def. 7.3 


0{s) 


ray direction in 7(5) 


Def. 7.3 


L+{s),L-{s),l+{s) 


positive, negative and effective ray length in 7(5) 


Def. 7.3 




paranietrization 


Def. 7.3 


Es 


boundary points with uniformly bounded ray lengths 


(7.37) 


a{s),f3{s) 


direction of ray and difference to tangent at 7(5) 


(7.44) 


m(s, •) 


mass coordinates 


(7.52) 


t(s,-) 


length coordinates 


(7.54) 


M(s) 


mass over 7(5) 


(7.53) 


E^,ei,ipi, 








corresponding quantities for a collection {71} 


Rem. 7.4 








EE,i,Os,i,'(pe,i, 








corresponding quantities for a collection {7e,i} 


Rem. 7.4 









2. Heuristics: the functional J^^ 

In the preceding section we postulated that the energy functional J^e, defined 
in (1.1), favours partial localization. More precisely, we claim that for small but 
fixed £ > functions {u, v) with moderate Ge{u, v) will be partially localized, in the 
sense that their support resembles one or more fattened curves. 

In this section wc first provide some heuristic arguments to support this claim, 
mostly to develop some intuitive understanding of the functional J^^ and its prop- 
erties. 

For the discussion below it is useful to remark that the distance di{u,v) as 
defined in Definition 3.1 scales as the mass J u = J v times a length scale, and can 
be interpreted as a u-weighted spatial distance between the mass distributions u 
and V. 

2.1. Fixed thickness (e — 1), part 1: comparing disc and strip structures. 

To gain some insight into the properties of J^^ we start with the case of fixed £ = 1, 
and study the limit of large mass 




We now compare two different specific realizations, a disc and a strip. 

(1) The disc: we concentrate all the mass of u into a disc in M^, of radius 
R M^/^, surrounded by the mass of v in an annulus (Figure 2). In this 
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{x : u{x) = 1} 
{x : v{x) = 1} 



1 2 1 



Figure 2. Disc and strip structures. 

setup mass is transported (in the sense of the definition of di{u, v j) over a 
distance 0{R) = 0(MV2), ^^^^ 

di{u, v)r~.MRr~. as M ^ oo. (2.1) 

On the other hand, the interfacial length / |Vu| equals 2'kR M^/^, and 
the functional J^i picks up the larger of the two: 

J^i{u,v) = di{u,v) + j |Vu| ~ 

(2) For the strip we take the support of u to be a rectangle of width 2 and 
length M/2, flanked by two strips of half this width for supp(t;) (Figure 2). 
The transport distance can now be taken constant and equal to 1, so that 
we find 

Ti{u,v) = M + M + 4 = 2M + 4. 

For this choice of geometry Ti scales linearly with M in the limit M oo. 

Comparing the two we observe that the linear structure has lower energy, for 
large mass, than a spherical one. This is a first indication that T\ may favour 
partial localization. 

With this example we can also illustrate an additional property. Let us take for 
supp(u) a strip of thickness t and length M/t, and recalculate the value of Tf. 

For large M this expression is dominated by the value of the prefactor t/2 + 2/t, 
suggesting two additional features: 

(1) There is a preferred thickness, which is that value of t for which t/2 + 2/t 
is minimal (i.e. t = 2); 

(2) For the preferred thickness, the functional Ti equals 2M plus 'other terms'. 

Although these statements only give suggestions, not proofs, for the case of general 
geometry, we shall see below that they both are true, and that these examples do 
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demonstrate the general behaviour. The 'other terms' in the example above are the 
single term 2t, a term which is associated with the ends of the strip; we now turn 
to a different case, in which there are no loose ends, but instead the curvature of 
the structure creates an energy penalty. 

2.2. Fixed thickness (e = 1), part 2: comparing line with ring structures. 

In the previous section we argued that the energy T\ favours objects of thickness 2 
and penalizes loose ends. We now consider ring structures, and we will see that 
also the curvature carries an energy penalty. 

Let the supports of u and v be the ring structures of Fig. 3: the support of u is 
a single ring between circles of radii ri and ra, and the support of v consists of two 
rings flanking supp(u). 



{x : u{x) = 1} 
{x : v{x) = 1} 




ri r2 



^3 ri 



\x\ 



ri r2 



rs ri 



\x\ 



Figure 3. Ring structures; the pairs {u,v) on the right-hand fig- 
ure indicate the values of u and v in that region in the plane. 



We first do some direct optimization to reduce the number of degrees of freedom 

in this geometry. For fixed outer radii ri and r^ one may optimize J^i over variations 
of r2 and ra that respect the mass constraint J u = J v, and this optimization results 
in the values of ri and rs given by 

r2 = i(i?2 + r?), and rl = l{R^+rl) (2.2) 

in terms of the mean radius R= ^ (ri +r4) . The functional !Fi can then be computed 

explicitly in terms of ri and r4 (see Appendix B). The interesting quantity is 
actually the energy per unit mass, J^i/M, as a function of the mean radius R and 
the thickness t of the structure, 

M:=Ju, 

Expanding TxjM around R = oo and t = 2,we find 

^iR,t) =2+\{t-2f + \R-' + Oi\t-2f + R-% (2.3) 
see Appendix B. 
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Again we recognize a preference for structures of thickness t = 2; we now also 
observe a penahzation of the curvature in the term i?~^/4. The main result of this 
paper indeed is to identify this curvature penalization for structures of arbitrary 
geometry, in the form of an elastica limit energy. 

2.3. Rescaling and renormalization. In the previous sections we have seen that 
the functional at £ = 1 has a preference for structures of thickness 2. It is easy to 
see, by repeating the arguments above, that this becomes a preference for thickness 
2e in the general case. For instance, the development (2.3) generalizes to 

J(i?, t)=2+ lit/e - 2f + \e'R-^ + 0{\t/e -2\' + e^R-'). 

This expression provides us with a recipe for identifying the bending energy in 
the limit £ — > 0. The curvature, which is equal to for this spherical geometry, 
enters the expression above as the term £^i?~^/4. This suggests that the alternative 
functional 

g,{u, v) := ^ [J^,{u, v) - 2M] , where M = ju = j v, 

may have a (Gamma-)limit similar to the Elastica functional. This suggestion is 
proved in Theorem 4.1. 

3. MONGE-KANTOROVICH distance, systems of curves AND A GENERALIZED 

CURVE BENDING ENERGY 

In this section we introduce some basic definitions and concepts. 

Definition 3.1. Consider u,v £ L^(M^) with compact support in satisfying the 
mass balance 

udC'^ = [ vdC^ = 1. (3.1) 

The Monge-Kantorovich distance di{u,v) is defined as 

di{u,v) := min( / \x - y \ d'j{x,y)) (3.2) 

where the minimum is taken over all Radon measures 7 on x with marginals 
uC^ and vC^, that means 7 satisfies 

I ip{x) d'y{x,y) = / ipudC"^, (3.3) 
/ ^P{y)dj{x,y) = [ i^vdC^ (3.4) 

jR2xR2 JR2 

for all C0(R2). 

The Monge-Kantorovich distance is characterized by the optimal mass transport 
problem described in section 5. 

Remark 3.2. For fj,,y &Vi, where 

Pi(R2) := Radon measure on R : / dfi — 1, \x\d^{x) < 00}, 

we can define the Monge-Kantorovich distance di{iJL,u) analogously to (3.2), sub- 
stituting uC"^ and vC"^ in (3.3), (3.4) by fi and u. The minimum in (3.2) is always 
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attained and di is a distance function on Vi [6]. Moreover di is continuous in the 
sense that for {iik)keN C Pi{M.'^), fj, e ri{R^) 



Hk ^ IJ- 3iS Radon measures and 
Jr2 \x\dnk /r2 \x\dfj,, 



see [6, Theorem 3.2]. 

A family of related distances is given by the p-Wasscrstein distances 

dp{ii,v) := min( / \x - yY" d'y{x,y)Y , (3.5) 

where the minimum is taken over all Radon measures 7 on x with marginals fj, 
and ly. The Monge-Kantorovich distance coincides with the 1-Wasserstcin distance. 

We will describe the limit structures in terms of systems of closed W"^'^ -curves. 
The following definitions are mainly taken from [12]. 

Definition 3.3. We call 7 : M ^ a closed W'^''^ -curve if 

7e<i(R,R^ 

7'(r) ^ for all r e R, 

7 is L-periodic for some < i < 00. 

// 7 is L-periodic and injective on [0, L) the curve 7 is called simple. The length 
1/(7) of a closed curve 7 is defined by 

L 



L(7):= / h'{r)\dr, 
Jo 



where L is the the minimal period ofj. 

A VK^'^-system of closed curves is a finite collection of closed W"^'"^ -curves. We 

represent a system of curves by a multiset {7i}i=i at, iV € N , i.e. a set in which 

repeated elements are counted with multiplicity. 

A system of curves T ~ {7;},-! n is called disjoint if each 7^ is simple and 

7j(R) n 7j(R) = for i ^ j. T has no transversal crossings if for any I < i, j < m, 
Si, Sj e R, 

7i(si) = 7j(sj) implies that jKsi) and^'j{sj) are parallel. 

We define the support, multiplicity and total mass of a system of curves T = 
{li}i=i,...,N by 

N 

supp(r) := IJ 7i(R), 



e{T,x) := #{{i,s) ■.ji{s)=x, l<i<N, 0<s<L{J^)}, 

N 

|r| 



and define a corresponding Radon measure nr on R^ to be the measure that satisfies 
/ ^dMr = Yl / y^{li{s))H{s)\ds (3.6) 

JR2 Jo 



for all tp€C[ 
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Two systems of curves are identified if the corresponding Radon measures coin- 
cide. 

Remark 3.4. We can represent a given system of closed curves F by a multiset 

r = {7i}i=i,...,Ar, where for all i — 1. .... 

7i is one-periodic, with 1 being the smallest possible period, (3.7) 
ji is parametrized proportional to arclength. (3.8) 

We generalize the classical curve bending energy defined in (1.5) to W^'^-systems 
of closed curves. 

Definition 3.5. Let T be a W'^''^ -system of closed curves represented by T = 
{7i}i=i,...,jv as in (3.7)-(3.8) and set Li = L{ji). Then we define 

W(r) := \Y,Lf j\'{{sfds. (3.9) 

Remark 3.6. The support and total mass of a system of curves coincide with the 
support and total mass of the corresponding Radon measure [46] . The multiplicity 
functions of F and of /ir coincide T^-'^-almost everywhere on the support of F. 
By Remark 3.8, Proposition 4.5, and Corollary 4.8 of [13], we obtain that F is a 
VF^'^-system of closed curves without transversal crossings if and only if /ip is a 
Hutchinson varifold with weak mean curvature H (E L^(/ir) such that in every 
point of supp(F) a unique tangent line exists. Moreover 

W(F) = \H\^diiT (3.10) 

holds. 

4. Main results 

To investigate the limit behaviour of the functionals .F^ as e — > we fix M > 
and study the rescaled functionals 

g,{u,v) := ^(T,{u,v)-2My (4.1) 

Theorem 4.1. The curve bending energy W as defined in (3.9) is the Gamma-limit 

of the functionals in the following sense. 

1. Let {u^,v^)^>o C L^(]R^) x L^{E?), R > and a Radon measure ji on he 
given with 

supp(?ie) C Br(0) foralleyO, (4.2) 

UgC^ ^ ji as Radon measures on M^, (4-3) 

and 

liminf C/e(Me, Wg) < oo. (4.4) 

£—►0 

Then there is a W^'^ -system of closed curves F such that 

2nr = 1^, (4.5) 

supp(F) is bounded, (4.6) 

F has no transversal crossings, (4.7) 

2|F| = M, (4.8) 
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and such that the liminf-estimate 

>V(r) < limmig.iue) (4.9) 

e— >0 

holds. 

2. Let T be a W'^''^ -system of closed curves such that (4.6) -(4.8) holds. Then there 
exists a sequence {us,Vs)e>o C tCg such that (4.2) holds for some R> 0, 

UgjC^ — > 2/ir as Radon measures on R^, (4-10) 

and such that the limsup-estimate 

W{T) > limsupae(we) (4.11) 

is satisfied. 

Remark 4.2. Let X be the space of nonnegative Radon measures on M^. Define 
a concept of convergence on X x X by which 



(/UijZ/i) {ii,v) mX X X iff 



Hi fj. as Radon measures on and 
U^>^^ supp(/ii) is bounded for some io € N . 

Note that no condition is placed on Vi. We now define 



and 



Qe{u,v) if = uC"^, V = vjC^ for some {u,v) G /Cg, 
oo otherwise 



>v(r) 



if H = v is given by a H^^'^-system of closed curves 
r = {7i}i=i,...,jv satisfying (4.5)-(4.8) 



oo otherwise. 
Then Theorem 4.1 can be rephrased as 

lim^e = ^0 (4.12) 

in the sense of Gamma- convergence with respect to the topology of X x X defined 
above. 

Remark 4.3. Implicit in the formulation of Theorem 4.1 is the following (contra- 
positive) statement: if a sequence {ue,Vs)e>o C ICg converges to a measure that 
cannot be represented by a system of curves satisfying (4.5)-(4.8), then is un- 
bounded along this sequence. Put differently, each of the following is penalized in 
at an order lower {i.e. larger) than 

(1) a fracture in the limit structure, occurring as a nonclosed curve; 

(2) a total mass |r| that deviates from M/2; 

(3) a non-even multiplicity. 

In the Introduction we conjectured the asymptotic development (1.6). Theo- 
rem 4.1 justifies the zeroth and second order of this development in the following 
way: if — 2M is of order 0(£^), then the three types of degeneracy above can 
not occur, and the dominant term is given in the limit by e^W. 
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Remark 4.4. By the definition of the set JCe and the functionals Te in (1.1), (1-2) 
it is clear that for a sequence (we,Ue)e>o along which C/e is bounded the supports 
of Ue and Ve necessarily concentrate as £ ^ 0. The crucial conclusions are (a) that 
the concentration is on curves rather than on points or other sets and (b) that the 
limit measure can be constructed as a sum of curves, each with density two. The 
latter property implies a uniform thickness of structures with bounded energy and 
paraphrases the stretching resistance of lipid bilayers. 

Remark 4.5. Wo also obtain a compactness result for sequences with bounded Q^'- 
Let {us,Ve)e>o C ^^(K^) X L^{M?) satisfy (4.2), (4.4). Then there exists a Radon 
measure on such that (4.3) holds for a subsequence e ^ 0. Theorem 4.1 then 
implies that there is a M^^'^-system F of closed curves such that the conclusions 
(4.5)-(4.9) hold. 

5. The Mass Transport Problem 

In this section we recall the classical Monge problem of optimal mass transport 
and review some results that we will use later. 

Definition 5.1. Fix two non-negative Borel functions u,v G L^(M.'^) with compact 
support satisfying the mass balance (3.1). By A{u,v) we denote the set of all Borel 
maps S -.M? ^R"^ pushing u forward to v, that is 

/ 'q{S{x))u{x)dx = / T}{y)v{y)dy (5.1) 
J Jm?- 

holds for all'neC°{R'^). 

Optimal mass transport. For u,v as in Definition 5.1 minimize the transport 

cost / : A{u, v) R, 

I{S) := / \S{x)-x\u{x)dx. (5.2) 

JQ 

The following dual formulation is due to Kantorovich [33]. 

Optimal Kantorovich potential. For u,v as in Definition 5.1 maximize K : 
Lipi(R2) ^ 

K{(t)) := [ (l){x){u - v){x)dx, (5.3) 

where Lipi(IR^) denotes the space of Lipschitz functions on with Lipschitz con- 
stant not larger than 1. 

There is a vast literature on the optimal mass transportation problem and an 
impressive number of applications, see for example [25, 48, 16, 5, 50, 32, 40]. We 
only list a few results which we will use later. 

Theorem 5.2 ([16], [27]). Let u,v be given as in Definition 5.1. 

1. There exists an optimM transport rn,ap S G A{u,v). 

2. There exists an optimal Kantorovich potential 4> G Lip2(M^). 

3. The identities 

di{u,v) = I{S) = K{ct>) 

hold. 
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4- Every optimal transport map S and every optimal Kantorovich potential 4> satisfy 
(j){x) — (j){S(x)) = \x — S{x)\ for almost all x G supp(u). (5.4) 

The optimal transport map and the optimal Kantorovich potential are in general 
not unique. We can choose S and ^ enjoying some additional properties. 

Proposition 5.3 ([16, 27]). There exists an optimal transport map S G A{u,v) 
and an optimal Kantorovich potential </> such that 

(j){x) = min {(j){y) + \x — y\) for any x e supp(u), (5.5) 
i/esupp(u) 

(^(y) = max (^(j){x) — \x — y'\) for any y G supp{v) , (5.6) 

xGsupp(u) 

and such that S is the unique monotone transport map in the sense of [27], 

+ ^rl^l^'^'i + for all x^^X2G with S{x^) ^ S^x^). 

\Xi-X2\ \S{xi) - S{X2)\ 

We will extensively use the fact that, by (5.4) the optimal transport is organized 
along transport rays which arc defined as follows. 

Definition 5.4 ([16]). Let u,v be as in Definition 5.1 and let (f) G Lip]^(]R^) be the 
optimal transport map as in Proposition 5.3. A transport ray is a line segment in 
M? with endpoints a,h G'M? such that (j) has slope one on that segment and a, b are 
maximal, that is 

a S supp(u), b G supp(ti), a b, 
(j){a) - (f){b) = \a-b\ 
\(l){a + t{a - b)) - (j){b)\ < \a + t{a -b)-b\ for all t > 0, 
\(f){b + t{b - a)) - (f){a)\ < \b + t{b - a) - a\ for all t > 0. 

We define the transport set T to consist of all points which lie in the (relative) 
interior of some transport ray and define £ to he the set of all endpoints of rays. 

Some important properties of transport rays are given in the next proposition. 

Proposition 5.5 ([16]). 

1. Two rays can only intersect in a common endpoint. 

2. The endpoints £ form a Borel set of measure zero. 

3. If z lies in the interior of a ray with endpoints a G supp(w),6 G supp(v) then (j) 
is differentiable in z with V(j){z) = {a — b)/\a — b\. 

In Section 7.1 we will use the transport rays to parametrize the support of u and 
to compute the Monge-Kantorovich distance between u and v. 

6. Overview of the proof of Theorem 4.1 

6.1. The lower bound. The heart of the proof of the lower bound (4.9) is a 

parametrization of supp(7i) U supp(u) that allows us to rewrite the functionals .?> 
and Ge in terms of the geometry of supp(m) and supp(t;) and the structure of the 
transport rays. 

The first step in the proof of (4.9) is to pass to functions u and v with smooth 
boundaries (Proposition (7.1)) so that i9supp(u) can be represented by a collection 
{7i} of smooth curves. For the discussion in this section we will pretend that there 
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is only one curve 7 : [0, L) — > R^, which is closed and which we parametrize by 
arclength. 

The properties of the Monge-Kantorovich distance di give that almost every 
point X e supp(u) U supp(w) lies in the interior of exactly one transport ray; the set 
supp(ii) U supp(u) can therefore be written, up to a null set, as the disjoint union 
of transport rays. 

In addition, along each ray mass is transported 'from u to v\ implying that each 
ray intersects 9supp(u). We therefore can parametrize the collection of rays by 
their intersections with 9supp(u). 

The ray direction of the ray TZ{s) that passes through 7(3) defines a unit length 
vector 9(s); we now introduce a parametrization ip by 



In terms of such a parametrization the first term of J> has a simple description as 
the length of the curve 7, 



To rewrite di [u, v) we use the fact that transport takes place along transport rays; 
when restricted to a single ray the transport problem becomes one-dimensional, 
and even explicitly solvable. 

The central estimate in this part of the proof is the following lower bound: 



In this inequality the geometry of supp(u) is characterized by functions a, (3, and M, 
which are illustrated in Figure 4. The functions a and [3 are angles: a{s) is the 



Figure 4. The parametrization of supp(u) U supp(ti). On the left 
six transport rays are drawn over the two supports; on the right 
the functions a, /3, and M are indicated. At a boundary point 
7(5), a{s) is the angle of the transport ray through that point, 
and (3{s) is the angle between 7 and the ray. Finally, M{s)ds is 
the amount of mass contained between the rays at positions s and 
s + ds. 



^{s,t) ■.^j{s)+t0{s). 





(6.1) 




d siipp u — J 
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angle between the ray direction 9{s) and a reference direction and /3(s) is the angle 
between 0{s) and the tangent 7'(s) to 7. The function M(-), finally, measures the 
amount of mass supported on a ray. The relation between the function M(-) and 
the scalar M in (1.2) is 

M = [ M{s)ds. 
Jo 

With the estimate (6.1), the main result can readily be appreciated. Let a 
sequence {ue,Vs) be such that 5e(tte,Ue) is bounded as £ ^ 0, and let a^, /?£, and 
Mg be the associated geometric quantities. With the inequality above, 

gs{us,vs) = 5 



> 







f 

Jo 



£2 £2sin/3e(s) A sin^ f]^{s] 

4sin5/3,(s) 

(6.2) 



ds 



Note that each of the three terms above is non-negative. If Qg (u, v) is bounded, 
then one concludes that as £ — > the mass M^{s) tends to one for almost all s. 
Similarly, converges to 7r/2, implying that the ray angle converges to the angle 
of the normal to 7^, and in a weak sense therefore a'^ converges to the curvature 
Ke = 7". The inequality above suggests a W^^^ bound for curves 7^ with the last 
term approximating 



This integral is the bending energy that we expect in the limit. 

Based on the inequality (6.2) we show that the boundary curves of supp(ue) arc 
compact, that the limit is given by a W"^'"^ system of closed curves as in (4.5)-(4.8) 
and that the lim-inf estimate is satisfied. We finally conclude that the limit of the 
mass distributions is identical to the limit of the boundary curves. 

The main reasons why the proof in Section 7 is more involved than this brief 
explanation are related to the gaps in the reasoning above: 

• A ray can intersect 9supp('u) multiple times, and care must be taken to 

ensure that the paramctrization is a bijcction. 

• The ray direction 9 is not necessarily smooth; to be precise, when the ray 
length tends to zero (which is equivalent to M(-) vanishing), 6 may vary 
wildly. 

• Similarly, when M(-) vanishes, the L^-estimate on a' in (6.2) degenerates, 
resulting in a compactness problem for the boundary curves. 

6.2. The upper bound. With the machinery of the lower bound in place, and 
with the insight that is provided, the upper bound becomes a relatively simple 

construction. Around a given limit curve tubular neighbourhoods arc constructed 
for the supports of u and v. The thickness of these neigbourhoods can be chosen 
just right, and the calculation that leads to (6.1) now gives an exact value for Qg. 
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The main remaining difficulty of the proof is the fact that a M^^'^-curve need not 
have a non-self-intersecting tubular neighbourhood of any thickness, and therefore 
an approximation argument is needed. 

7. Proof of the lim-inf estimate 

In this section wc prove the first part of Theorem 4.1. The main idea is to 
use the optimal Kantorovich potential of the mass transport problem of Ug to Vg 
to construct a parametrization of the support of and and to derive a lower 
bound for the functionals 0^{u^,v^). This lower bound yields the compactness of 
the systems of curves which describe the boundary of the support of Ug as well as 
the desired lim-inf estimate. 

We first show that we can restrict ourselves to a class of 'generic data'. 

Proposition 7.1. It is sufficient to prove the liminf part of Theorem 4-1 under the 
additional assumptions that 

M=l, thatis ju. = jv. = l /orK,..)./C., (7.1) 
su'pQe{ue,Ve) < A for some < A < oo, (7.2) 

e>0 

{ue,Ve) e /Ce for all £ > 0, (7.3) 
YmiQ^{u^,v^) = liminf ^£ («£,«£ ) exists, (7.4) 

e— >0 £— >0 

9supp(u£) is given by a disjoint system of smooth, simple curves F^. (7.5) 

Proof. To prove that we can restrict ourselves to M = 1 we do a spatial rescaling 
and set 



e := Ue{x) := Mu^{Mx), Vg{x) := Mv^iMx). 



We then obtain 



I'- /./' = '• .-/lVS.I^^/|V«.|. (7.6) 

Moreover, if S is the optimal transport map from to v^, then S{x) := M^^S{Mx) 
is the optimal transport map from ug to vg and 

dl{us,Ve) = M-'^di{Ue,Ve). (7.7) 

We obtain from (7.6), (7.7) that 

= Mg,iu„v,). (7.8) 

We observe that {ug, vg) G JCg iff {us,Ve) e tCg and that {ug)g>o satisfies (4.2), (4.3) 
for R := R/M and // defined by 

J v{x) df,{x) = j ^V{^) dn{x) for v € C°{R'). (7.9) 
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If Theorem 4.1 holds for M = 1 then there exists a M^^'^-system of closed curves F 
with 

= fi, 2|f| = 1, (7.10) 

such that (4.6), (4.7) holds and 

W(f) < limMgg{u,,v,). (7.11) 

e— »0 

We deduce from (7.9) and (7.10) that /x is given as system of closed curves: Let 
f = {7i}i=i,...,jv with 

and define a system of curves F = {7i}i=i,...,Ar by 

-fi:[0,Li) ^ R^, L,:=ML„ 
■■= forse[0,L,). 
Then |7,'| = 1 and F satisfies (4.6), (4.7). We obtain from (7.9) that 

/p ^" pLi 

rjdfi = M J ri{Mx)djl{x) = 2A/^ j r]{M^i{s))ds 

m -Li r 

= 2J2 ■n{^i{s))ds = 2 -ndnT, 

which yields ji = 2/xr and, by taking ri = 1 inside -Bi{(0), that (4.8) holds. Finally 
we deduce from (7.9) 

^ i=l -^0 

and by (7.8), (7.11) this yields that (4.9) holds. Therefore we have shown that it 
is sufficient to prove the lim-inf part of Theorem 4.1 for M = 1. 

Eventually restricting ourselves to a subsequence £ we can assume (7.4). 
The existence of a < A < oo such that (7.2) holds follows from (4.4). By the 
definition of Ge this in particular implies (7.3). 

It remains to show that it is sufficient to consider functions such that 
3supp(u£) is smooth. We already have seen that we can assume that {ug, Vg) € ICg. 
This implies that supp(u£) is a set of finite perimeter. Moreover supp(u£) C Bii{0) 
and by [7, Theorem 3.42] there exists a sequence of open sets {Eh)keN with smooth 
boundary such that 

Xe, eue inL\R^), (7.12) 

^ e|Vu,|(R2) (7.13) 
Ek C B2r{Q). (7.14) 
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We set 



and define functions Uk, 



Tk := e^\Ek\-^ (7.15) 



Uk{x) := IxE^rk'^)- (7-16) 



This yields 



and 



/ 



e 



Uk = ^-^\Ek\ = 1 (7.17) 



[ \Vuk\ = - f |V£fe|. (7.18) 

Moreover, from (7.12), (7.15) and | supp(w£)| = £ we observe tliat 

rfe ^ 1 as fc — > oo. (^-IQ) 

We compute that 



/ \uk-Ue\ = - [ ^Ek{rk^3:) - SUe{x) dx 

< 7 / \^Eu{rk^x)-eUe{r^^x)\dx 



+ - / \eue{r^.^ x) — eue(x)\dx 



£ ./b2 

„2 



~ / \XEk - sus{x)\dx + / \ue{r^^x) - Us{x)\dx. 



Since the right-hand side converges to zero as fc — > 00 by (7.12), (7.19) we deduce 
that 

Uk -> We mL^{R^) (7.20) 
as fc ^ 00. Equations (7.13), (7.18) and (7.19) yield that 

/ IVufel / |Vue| asfc^oo. (7.21) 

We construct approximations Vk of Vg by setting 

Vk{x) := Vs{x){l-suk{x)) +Wk{x), (7.22) 
for a suitable Wk S L^[B2r{0), {0, £~^}) with 

{wk - evsiik) = 0. (7.23) 



We then deduce from (7.22), (7.23) that {uk,Vk) G ICe- It follows from (7.20), (7.23) 
and UeVe = that ^ as A: ^ 00. Together with (7.22) we obtain that 

lim Vk = Ve in ^^(R^). (7.24) 

k — *oo 

By (7.14), (7.20) and (7.24) we deduce from Remark 3.2 the continuity of the 
Monge-Kantorovich distance term in Ts{uk,Vk)- Together with (7.21) this yields 
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Therefore we can choose k{e) G N such that for Ue :— Uk{e)t '■= Vkie) satisfy the 
estimate 

\\Ue-Ue\\L^R2^ + \Qe{Ue,Ve)-Gs{Us,Ve)\ < £■ 

This yields a sequence {ug,Vs)s>(i that satisfies (7.5) and has the same limit as 
(we,v.),>o in (4.3), (4.9). □ 

Prom now on, for the rest of this section, we assume that (7.1)-(7.5) hold. 

7.1. Parametrization by rays. Due to (7.5) the boundary 9supp(u£) is smooth. 

We arc free to choose a suitable representation. 

Remark 7.2. We can represent c?supp(M£) by a disjoint system of closed curves 
such that for all 7 e 

7 is parametrized by arclength, |7'| = 1, 

supp(u£) is 'on the left hand side' of 7. (7-25) 

In particular 7 is Z/(7)-periodic and 

det(7'(s),z/(s)) < 0, 
where y{s) denotes the outward unit normal of supp(u£) at 7(5). 

Let (j) e Lipj(lR^) be an optimal Kantorovich potential for the mass transport 

from ti^ to as in Proposition 5.3, with T being the set of transport rays as in 
Definition 5.4. Recall that (j> is differentiable, with |V</)| = 1, in the relative interior 
of any ray. 

Definition 7.3 (Parametrization by rays). For 7 € we define 

(1) a set E of 'inner points ' with respect to the transport set T , 

E := {s e K : 7(5) e T}, 

(2) a direction field 

e:E^S\ 6{s):=V4>{l{s)), 

(3) the positive and negative total ray length L+, L~ : £ — > M, 

L+{s) := sup{t > : ^(7(5) + te{s)) - </.(7(s)) = t}, (7.26) 

L-{s) := m{{t < : .^(7(5) + t6{s)) - <j>{^{s)) = t}, (7.27) 

(4.) the effective positive ray length : M, 

Z+(s) := sup {t e [0, L+{s)] : 7(3) + t6I(s) e Int(supp(ue)) for allO<T<t}. 

(7.28) 

Finally we define a map ip which will serve as a parametrization of supp(M£) U 
supp(w£) hy 

ip : [0,L(7)) X M ^ IR^ V(s,*) := 7(s) +i6'(s). (7.29) 

Remark 7.4. All objects defined above are properties of 7 even if we do not 

denote this dependence explicitely. When dealing with a collection of curves {7^ : 
i = I, ■■■,N} or {7e^j : e > Q,i = 1, ...,N{e)} then Ei, l'^^ etc. refer to the objects 
defined for the corresponding curves. 
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Figure 5. The parametrization in a 'simple' situation. Here 

l~^{s) = L^{s) coincide. 

See the Figures 5 and 6 for an iUustration of the parametrization defined above. 

The effective ray length is introduced to obtain the injectivity of the parametrization 
in the case that a ray crosses several times the boundary of supp(M£) (as it is the 
case in the situation depicted in Figure 6). The set {7(5) : l~^{s) > 0} represents 
the points of the boundary where mass is transported in the 'right direction'. 
We will first prove some results which are analogous to those in [16]. 

Lemma 7.5. The ray direction 9 : E ^ M.'^ is continuous. The positive, negative, 
and effective positive ray lengths L"*", L~ , : £J — » R are measurable. 

Proof. Consider a sequence (sfe)fegN C E, s ^ E, with Sk s. Choose tk < tk 
such that 

-'^<tk < L+{sk), (7.30) 

l+{sk) - i < 4 < l+{sk). (7.31) 
It follows from (7.26) and (7.30) that 

<f>{l{sk)+tke{sk)) = cj,{j{sk))+tk. (7.32) 

Since {tk)keN , (ife)feeN are uniformly bounded by 2R and |^?(sfe)| = 1 there exists 
a subsequence k ^ 00 and f, i e M, (9 e with t < t, |^^| = 1 such that 

tk ^ t, ik ^ i, and 0{sk) 9. 

We deduce from (7.30)-(7.32) that 

lim sup L+(sfe) < t, (7.33) 

lim sup Z+(sfc) < i, (7.34) 

fe— s-OO 

<t>{l{s)+t9) = <l>{i{s))+t. (7.35) 
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Since <j) has Lipschitz constant one it follows from (7.35) that 

e = V(/)(7(s)) = e{s). (7.36) 
We therefore deduce by (7.35) that t < L^{s) and by (7.33) we find that 

limsupL~''(sfc) < L~^{s), 

k—^oo 

which proves the upper-semicontinuity and therefore the mcasurability of . By 
analogous arguments one proves that L~ is lower-semicontinuous and measurable. 
Let now < r < f . Since ik ^ i we obtain from (7.31) that 

4'{l{^ki) + TO{ski)) S supp(u£) for all sufficiently large A; e N . 

Since supp(u£) is closed and since <!) is continuous this yields 

(pl^is) + t9{s)) e supp(u). 

But < T < f was arbitrary and we deduce that i < l~^{s). By (7.34) this shows 
the upper-semicontinuity and measurability of Z"*". □ 



□ {.x : n{x) = 1} 

□ {x : vix) = 1} 




Figure 6. A situation where a ray crosses several times the 
boundary of supp(we), given by curves 71,72,73- On the right 
for each of the curves 7j the positive, negative and effective ray 
lengths L^,i^,Z^ is depicted. Observe that ^2^(5) = 0. 

The next result is similar to Lemma 16 of [16]: since rays can only intersect in 
a common endpoint, the ray directions vary Lipschitz continuously. 

Proposition 7.6. Let 7 € fee given as in Remark 7.2 and set for 5 > 

Es := {sGE: L+{s), \L-{s)\ > S}. (7.37) 
Then 6 is Lipschitz continuous on Es with 

\e{s,) - e{s2)\ < 7IS1-S2I for allsx,S2€E5. (7.38) 

Proof. To prove (7.38) it is sufficient to consider the case 0{s\) ^ 9{s2). Then there 
exist t\,t2 € M such that 

7(si)+ii^(si) = 7(s2)+i2^(s2). (7.39) 
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By the definition of L"*", L~ and since rays can only intersect in common endpoints 
we obtain 

ti &^\{L-{si),L+{si)), i = l,2. (7.40) 

By taking the determinant with 6(82) in (7.39) we compute 

det (7(si) - 7(S2), e{s2)) = -ti det (^(si), ^(sa)) . (7.41) 

Prom (7.40) and since si, S2 € Eg we deduce that > S^. Taking squares in (7.41) 
we therefore obtain 

(S1-S2)' > 5'det(^(si),e(s2))' 



(52 



i-{e{s,).e{s2)f 

e{s,)-e{s2)f\9{s,) + e{s2)\". (7.42) 



Next we observe that 

\0{si) + e{s2)\ > 1 if |si-S2| <(5. (7.43) 

Indeed we compute 

4^ = </.(7(si) + S9{s,)) - 0(7(.si) - Se{s,)) + ^(j{s2) + 66(32)) - Hl{s2) - (50(s2)) 
= (/.(7(si) + 5e{si)) - 0(7(52) - 6e{s2)) + <^(7(S2) + 5e{s2)) - <t>{l{si) - 69{si)) 

< 2|7(si)-7(s2)|+2(5|^?(si)+0(s2)| 

< 2\si ~ S2\ + 2S\e{si) + e{s2)\ 

which proves (7.43). 

By (7.42), (7.43) we deduce that (7.38) holds for si, S2 e Es with \si-S2\<S. In 
the case |si — S2\ > 6 we obtain (7.38) since \9{si) — d{s2)\ < 2. □ 

Definition 7.7. Define two functions a,p : E — > R mod27r by requiring that 

0{s) = (^^2[l)) ' det(7'(s),0(s)) = sin/3(s). (7.44) 
We note that 

sin/3(s) > hi seE with l+{s) > (7.45) 

by (7.25) and the definition of Z+. 

Proposition 7.6 yields the Lipschitz continuity of 9, a on sets where the posi- 
tive and negative ray lengths are strictly positive. By Kirszbraun's Theorem and 
Rademachcr's Theorem wc obtain an extension which is diffcrcntiablc almost ev- 
erywhere. To identify these derivatives as a property of a we prove the existence 
of a' in the sense of approximate derivatives [26, section 6.1.3]. 

Lemma 7.8. The function a : £ — > M mod 27r as defined in Definition 7. 7 is almost 
everywhere approximately differentiable and its approximate differential satisfies 

sin f3{s) -ta'{s) > for all L-{s) <t < L+{s) (7.46) 

for almost all s G E with l~^{s) > 0. 
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Proof. Consider for (5 > the set Es defined in (7.37). By Proposition 7.6 the 
restriction of 6 to Eg is Lipschitz continuous. From Kirszbraun's Theorem we 
deduce the existence of a Lipschitz continuous extension ci^ : R ^ M mod 27r of 
a, see for example the construction in [19]. This extension is almost everywhere 
differentiable by Rademacher's Theorem. By [26, Corollary 1.7.3] almost all points 
of Eg have density one in Eg and we deduce that a is approximately differentiable 
in these points and that the approximate differential coincides with a'g. Since 
E = UfceN -E'l/fe this proves that a is approximately differentiable almost everywhere 
in E. 

Since a measurable function is almost everywhere approximately continuous by 
[26, Theorem 1.7.3], it is sufficient to prove (7.46) for s e £ with l+{s) > such 
that 

• a is approximately differentiable in s, and 

• there exists a sequence {sk)keN C E with — > s as fc — > oo and 

L+{sk)^ L+{s), L-{sk)^L-{s), l+{sk)^l+{s), 
L+{sk), -L-{sk), /+(sfc) > for all fc e N . 

Since (7.46) is satisfied for a'{s) = we can also assume that 

• a'{s) > (the other case is analogous); 

• {a{sk) - a{s))/{sk - s) ^ a'{s), or equivalently, that {9{sk) - 9{s))/{sk - 
s) ^ e'{s); 

• for all A; e N there exist tk, t^. such that 

7{sk) + tke{sk) = 7(s) + tms). (7.47) 
Analogously to (7.40), (7.41) we find from (7.47) 

tk G M\ (i-(,sfe),i+(,sfe)), tl e R\{L-{s),L+{s)), (7.48) 

det (7(sfc) - 7(s), 0{s)) = -tk det {e{sk), 0{s)) . (7.49) 
Writing equation (7.49) as 

det(2M^,,(,)) = -i,det(^if^l^,.(.)), 

\ Sk — S / \ Sk — s / 

we observe that the left-hand side converges to det(7'(s), 6'(s)) = sin/3(s) > 
and the determinant on the right-hand side to dct(6' {s),9{s)) = —a'(s) < 0. For 
sufficiently large k we can therefore eliminate the negative range for tk in (7.48) 
and obtain 

tk > L'^isk)- 

Again using the negative sign of det((^(sfe) — 6{s))/{sk — s),6{s)) we then find 

det(2M^,^(.)) > -L+(..)det(^if^l^,^(.)) 

Vsfe — s / \ Sk — s J 

Taking the limit k oo we deduce 

sin/3(s) > L+(s)a'(s). 
This proves the Lemma. □ 

For each 7 e we define analogously to (7.29) a map ^. Restricting these maps 
suitably we obtain a parametrization of supp(u£) which is essentially one-to-one. 
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Proposition 7.9. Let = {'^i i — 1, ...,-/V} he chosen as in Remark 7.2. For all 
i = l, N let Li = i(7i). Then, with 

Di := {{s,t):seEin[0,Li),0<t<lt{s)}, (7.50) 

the restrictions ipi : ^ give, up to a Lebesgue nullset, an injective map onto 
supp{ue): for almost all x £ supp(u£) there exists a unique i £ {1,...,A''} and a 
unique {s,t) € Di such that ipi{s,t) = x. 

Proof. For both the injectivity and the surjectivity it is sufficient to consider only 
interior points of supp(ue) that are not ray ends, since the boundary of supp('U£) 
and the sets of ray ends form an £^-nuUset (Proposition 5.5). 

For the surjectivity, let x be such a point in supp(u£). By (5.5) there exists 
y e supp(v) such that 

(j){x) - (p{y) = \x-y\ 

and it follows that x is on the interior of a ray with direction W(l){x) = jfE^ • Define 

t := min {t > : x — tW4>{x) G dsupp{us)}. 

Then there exists i G {1, N} and s G Ei such that l~^{s) >t>0 and 

X = "fi{s) + tOi{s) = il)i{s,t). 

To prove the injectivity part, let x be the same point again, and assume that 
there exists j G {1, ...,iV}, (a, r) G Dj such that 

x = ^i{s,t) = 'ipj{a,T). 

Since different rays cannot intersect in interior points such as x, the three points x, 
7i(s), and 7j(cr) are on the same ray, and since t > and r > we have 

(t>i'yi{s)) < <i>{x) and <?i(7j(o-)) < <i>{x)- 

This implies that 7i(.s) and 7j((T) arc on the same side of the ray with respect to x. 
On the other hand by the definition of l^ and since t < l^{s), t < l^{(7), neither 
the part of the ray between x and 7i(s) nor the part between x and contains 
a point of 9supp(?y,j). This implies that 7i(s) = 7j(o') which proves the injectivity 
part of the proposition. □ 

Using Proposition 7.6 and Proposition 7.9 we can justify the following transfor- 
mation formula for integrals. 

Lemma 7.10. Let = {7i}i=i,...,jv be as in Remark 7.2 and Di as defined in 
(7.50). Then for all g G L\R'^) 

r ^ r 1 

g{x)ue{x)dx g{ipi{s,t)) - (sin Pi{s) -ta'iis)) dtds (7.51) 

holds. 

Proof. We deduce from Proposition 7.6 that for alH G {1, ...,iV} the functions 6i 
and ipi are approximately difFcrcntiablc on Di with approximate differentials 

' «\ ^ \^ COSQ!i(s) J 

d,iPi{s,t) = ^'i{s) + te'i{s), dtMs,t) = e'iis). 
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This yields that 

\det DMs,t)\ = I det +tdet {e'i{s),0,is))\ 

= |sin/3i(s) - ia-(s)| 

= sin/3i(s) - ta'iis) 
where we have used (7.46) in the last equality. 

We then deduce from the generalized transformation formula [8, Remark 5.5.2] that 
/ g{tpi{s,t))- (sin l3i{s) — ta^{s)) dtds = ( g{x)ue{x) dx. 

Summing these equalities over i = 1, N we deduce by Proposition 7.9 that (7.51) 
holds. □ 

Often it is more convenient to work not in length coordinates but rather in mass 

coordinates which arc defined as follows. 

Definition 7.11. For 7 G Fg and s G R we define a map trig : M ^ M and a map 
M :R^R by 

m(M) := ^/^"(^)>0' (7.52) 

1 otherwise. 

M{s) := m{s,l+{s)). (7.53) 



9supp« = 7 



Figure 7. The function m. Two rays delimit a section of w-mass 
(below the curve 7, in this figure) that is transported to the asso- 
ciated section of w-mass. For given t > 0, m{s, t) cLs is the amount 
of w-mass contained between the two rays and stretching from 7(5) 
to the point 7(s) + t6{s). For t < 0, m.{t, s) is the corresponding 
amount of t;-mass, with a minus sign. 



Lemma 7.12. Themapm{s, 
i{s,m) 



■ strictly monotone on (L (s), L^{s)) with inverse 



sin/3(s) 
a'{s) 



2a'{s)e 
sin^ /3(s) 



rnj ^ . 



(7.54) 
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Let Tg = {7i}i=i,...,w be given as in Remark 7.2, set Li = L{'ji), and let Mi be the 
mass function (7.53) for curve 7,. We obtain for g G that 

f ^ pLi pMi{s) 

J g{x)ue{x)dx = 'Y^J J g{tpi{s,ti{s,m)))dmds. (7.55) 

In particular, the total mass of is given by 

/ Ue{x)dx = Y] I Mi{s)ds. (7.56) 

Proof The calculation of the inverse is straightforward. For the transformation 
formula (7.55) we combine (7.51) with the remark that 
d 1 

— m(s,i) = -(sin/3i(s) -ta'^{s)). 

Note that the integration interval < t < /^(s) in the definition (7.50) of Dj 

transforms to < m < Mi{s). □ 

Typically, the optimal transport map is built by gluing solutions of one- 
dimensional transport problems on single rays together. We observe a similar 
structure here. 

Proposition 7.13. Let 7 €E Fg be as in Remark 7.2 and S be the optimal transport 
map from Ug to v^, as in Proposition 5.3. Define for s € E with l^{s) > an 
interval I{s) C M and measures f^,f^ on I{s), 

I{s) := (m,(L-(s)),m,,(L+(.))), 

df^ := Ue{'il^{s,t{s,m))) dm, df~ :— Ve{i^{s,t{s,m))) dm. 
Moreover, define a map 

5:/(s)nsupp(/+) ^ /(s)nsupp(/7), 

by the equation 

5(V(s,t(s,m))) = ^(s,t(s,5(m))) /or m e / n supp(/+). (7.57) 

Then for almost all s & E with l~^{s) > the map S is the unique monotone 
transport map pushing /+ forward to f~ . 

Proof. It is sufficient to prove the claim for almost all s € {l~^ > k~^}, k £ N 
arbitrary. Let A C E he any set such that 

Ac {sgE: l+{s) > k-^}, diam(A) < fc"^ (7.58) 

We define rays and a set of rays 

n{s) := {7(s)+T^(s):re(L-(s),i+(s))}, 

n{A) := y n{s). 

seA 

Since different rays can only intersect in a common endpoint, two rays TZ{s),TZ{a) 
with s,a G A, s ^ a can only be disjoint or identical. The latter case is excluded 
by (7.58): since (s) , l'^ {a) > kr^ the distance between 7(5) and 7((7) has to be 
at least fc^^ in contradiction to the assumption on the diameter of A in (7.58). 
Therefore the union {TZ{s) : s e A} is pairwise disjoint. 
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Let now (3j)jGN be a dense subset of C°(Ili). Since S pushes UsItz^a) forward to 
Ve\Ti{A) we obtain for all j € N 

/ gj{(l){S(x)))ue{x)dx = I gj{(p{x))vs{x)dx. (7.59) 
Jtz{A) Jti(A) 

Repeating the changes of variables of Lemmas 7.10 and 7.12, and using the fact 
that {7?.(s) : s € A} is pairwise disjoint we obtain 

^ I I { ^^^i'^^^^^'^^^''^^^^)^^^^^^'^^^'^^^^ dmds, 

and by using the definition of S (7.57) and the linearity of ^ along rays we write 
this as 



/ / gj{(t>{'y{s)) + t{s,S{m)))us{ijj{s,t{s,m))) dmds 

J A Jl{s) 

= 5j (0(7(5)) +t(s>w))we(-0(s,t(s,m))) 

J A Jl(s) 



I{s) 

Since A as above was arbitrary we deduce that for almost all s £ {Z+ > 0} and all 

/ gj{cl>{l{s))+t{s,S{m)))df+{m) = [ gj{^{^{s)) + t{s,m)) df- {m). 

Jl(s) Jl(s) 

Since {gj)j^^ is dense in C°(R) and since for fixed s the function m 1— > ^{'y{s)) + 
t{s, m) is a homeomorphism it follows that 

/ giSim))df+{m) = [ g{m)df-{m) 

Jl(s) Jl(s) 

for all g G C°(M). Therefore S pushes /+ forward to /~. The monotonicity of S 
follows from the monotonicity of S. By [5, Theorem 3.1] the map S is the unique 
monotone transport map from /+ to f~. □ 

Lemma 7.14. Let 7 e as in Remark 7.2. Then for almost all s G E with 
l~^{s) > we obtain that 

\^{s,t{s,m)) - S{il>{s,i{s,m)))\ > t{s,m) - i{s,m - M{s)) (7.60) 

for alio <m< M{s). 

Proof. We let S be as in (7.57). Using the monotonicity of t(s, .) we obtain 
|V'(s,t(s,m)) - S{tp{s,i{s,m)))\ = \tp{s,i{s,m)) - il){s,t{s,S{m)))\ 

= t{s,m)-t{s,S{m)). (7.61) 

By Proposition 7.13 the map S is the unique monotone transport map from /+ to 
/~. In particular 

f-{[S{m),S{M{s))]) = f+{[m,M{s)]) (7.62) 
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holds. Since restricted to the set [0, M{s)] coincides with the Lebesgue measure 
and since < /g", /+ < £^ we deduce from (7.62) that 

S{M{s)) - S{m) > M{s) - m (7.63) 

for any < m < M(s). The value S{M{s)) is less or equal than zero since 
S is decreasing and maps onto the support of f^. This implies by (7.63) that 
S{m) < m — M{s). By the montonicity of tg wc deduce 

t{s,m) -t{s,S{m)) > t{s,m) - i{s,m- M{s)). 

Together with (7.61) this proves (7.60). □ 

By Lemma 7.6 the ray directions vary Lipschitz continuously on sets where the 
positive and negative ray lengths are bounded from below. We obtain a Lipschitz 
bound also on sets where M(-) is bounded from below. 

Lemma 7.15. Let 7 e be as in Remark 7.2. Then the inequality 

sin/3(s)min(i+(s),|L"(s)|) > |m(s) (7.64) 

holds. In particular, for all < k < 1 the function a is Lipschitz continuous on the 
set {s : M{s) > (1 - k)} and 

lo'WI < («5) 
holds for almost all s gR with M{s) > (1 — k). 

Proof. By (7.46) we get that sin/3(s) > for s e with l+{s) > and that 
a'{s) = if sin/3(s) = 0. Since in the latter case (7.64) holds it is sufficient to 
assume that sin/3(s) > 0. We consider first the case a'{s) < 0. Prom (7.52) and 
(7.53) we obtain 

M{s) = -(^l+{s)smp{s)-^l+{sfa'{sj^ 

By Proposition 7.13 the estimate 

M{s) = /+((0,M(s)) < f-{l{s)n{L-is),0]) 
holds and we deduce that 



.0 -, 

M{s) < / - {sin j3{s)-ta'{s))dt 



iL-(s) 

< L^lMsin/J(s) 

e 

= (7.67) 
Together with (7.66) we deduce 

Therefore ^ = 212^M/+(s) satisfies the inequality 

f + 2M(s)C - 2M{sf > 
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and we obtain that 

> (V3-l)M(s). (7.68) 

e 

Together with (7.67) and l'^(s) < L^(s) this proves (7.64) in the case that a'{s) < 0. 
Next we assume that a'{s) > 0. From (7.52), (7.53) we observe that 

M{s) < ^L+{s) sin l3{s). (7.69) 

To prove (7.64) also for L~{s) let us assume — without loss of generality — that 
^ sin^(s) ^('^)- We then deduce that 

\L-{s)\ < L+{s) (7.70) 

and compute 

M{s) < / Ve{i){s,t)){sml3{s)-ta'{s))dt 

Jl-(s) 

< f - {sin P{s) - L- {s)a' {s)) dt = ^-^^-^ {sin /3{s)-L-{s)a'{s)) . 

Jl-{s) ^ ^ 

By (7.70) and (7.46) we deduce that 

M{s) < i^llM(sin/3(5)+L+(s)a'(s)) < 2|L-(s)| 

e e 

Together with (7.69) this proves (7.64) in the case a'{s) > 0. 

The estimate (7.65) follows from (7.64) and (7.46). □ 

7.2. Estimate from below. Let S be the optimal transport map and </> € 
Lipj^(M^) be an optimal Kantorovich potential for the mass transport from Us to Ve 
as in Proposition 5.3. 

Using Lemma 7.12 and Lemma 7.14 we obtain for the Monge-Kantorovich dis- 
tance between u to u that 



§io Jo 

N ,.L, fAh{s) 



di{u,v) = / \x — S{x)\us{x) ( 

\ipi{s,ti{s,m)) - S{^|Ji{s,ti{s,Tn))\ dmds 

^ z2 / [U{s,m)-U{s,m-Mi{s))]dm. (7.71) 
i=l -^0 ^0 

We further estimate the right-hand side of this inequality. 

Lemma 7.16. Let Tg be a disjoint system of curves parametrized by arclength as 
in Remark 7.2 and let 7 e F^. Then we obtain for all s G E with l~^{s) > and 
sin/3(s) > that 

f.M(s) 



/ [i{s,m) -i{s,m- M{s))]dm 
Jo 



^ -M(s)2 + — 4-— a'(s)2M(s)4 + i?(s)£5, (7.72) 



sin/J(s) ' ' 4sin5/3(s) 
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where 



< R[s) < l^^a'isr. 
- ^ ^ - 9sin/3(s)9 ^ ^ 



(7.73) 



Proof. We compute for r e 



t{s,rn)dm = 



sin/3(s) 





V sin^ Bis) / 



siii^/3(s)r/ 2a'{s)e \ 
3a'(s)2£ [V ~ sin^ /Jfs) / 



+ n r - 1 



sin^/3(s) / sin^/3(s) 



(7.74) 



and thus 



.Mis) 

/ [t(s,m) - t(s,m- M(s))] dm 
Jo 

= / i{s,rn)dm + 
Jo Jo 



t(.s, rn) dm 



sin^ /3(s) 
3a'(s)2e 



/ 2a'(s)e ,\i 2a'(s)e aI 

(l+ . 2 J , ^(^)) . 2 J , ^(^)) -2 

V sm li s) / V sm Ills) / 



(7.75) 



sin^/3(s) ' V V sin''/3(s) 
By a Taylor expansion wo observe that for all ^ > 

(i+e)^ + (i-o^-2 = le + ^^e+l'^ 

with < C(0 < ^- Using this expansion in (7.75) we obtain (7.72), (7.73). □ 



(7.76) 



Putting all information together we derive the following estimate. 

Proposition 7.17. Lei 9supp(u£) be given by a disjoint system of closed curves 
parametrized by arclength as in Remark 7.2. Then the lower bound 



Q.K.v.) > g£'[^(_L_-,)M.w^4{M.W-i) 



ds 



^io 4smA(s) Vsm/3i(s)/ 



holds. 



Proof. Since L, = L(7i) and by (7.56) we obtain that 



JV 



; / |Vue| = and 1 = / = / M,{s)ds. 

Jm Jo 



i=l ' 



Prom (7.71) and Lemma 7.16 we therefore obtain 



+E/ t^-^Vtt^^^^)'^'^^^'' 

4sin°A(s) 
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and observing that 

the estimate (7.77) follows. □ 

7.3. Compactness of the boundary curves and proof of the lim-inf esti- 
mate. In this subsection we prove that the boundaries of supp(u£) converge to a 
system of closed curvcis which satisfies the lim-inf estimate. We later identify this 
limit with the limit of u^C^ ■ 

Let 9supp(Me) be given by a system of closed curves = {pie,i\i=\,...,N{e) 
parametrized by arclength as in Remark 7.2. In particular we recall that 

The total length of the boundary of supp(u£) is given by IF^] and we obtain from 
(7.2) and (7.77) 

ir^l = V ie.i = £ / IVuel < 2 + e25,(we,^^e) < 2 + e2A. (7.78) 

This implies the convergence of the measures e|Vue| = [iVe Radon measures on 
for a subsequence e ^ 0. 

We first define modified curves 7^^, which are uniformly bounded in W^'^ (O, Lg^j) . 

Definition 7.18. Let 7 G Fe and < A < 1. Choose a periodic Lipschitz continu- 
ous function a^^i : M. mod 27r with period L^,i such that 



ae,i = a^,i on {s € M : M^,i{s) > 1 - A}, (7.79) 
2 



and set 



ysin a^^i J ' \^— cosa^,, 

We then define ^^^,1 : K — > to he the curve which satisfies 

leM = 7...(0), (7.81) 
I'eM = for all sew. (7.82) 

Remark 7.19. By Proposition 7.6 the function ae,i is Lipschitz continuous on the 
set {Me,i > 1 — A}. By Kirszbraun's Theorem and (7.65) an extension cts^i as in 
Definition 7.18 exists. 

The curves 7£,j are not necessarily closed (or, equivalently, the functions ^e,i are 
not necessarily periodic); the restriction ^e,i\[o,Le,i) may also have self-intersections. 
However the tangents of these curves, the functions 7^,, are periodic with pe- 
riod Le,i. 

Lemma 7.20. Let £>0, 0<A<1, l<p<2, and consider for i e {1, ...,N{e)} 
functions ctg^i and modified curves 7e,, as in Definition 7.18. Then there exists a 
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(7.83) 
(7.84) 



constant C(A, A), which is independent of s, such that 

E/ KM-t,i{^)\ds< eC{X,A) 
i=i Jo 

E / \l'JM'ds<C{X,A), 

i=l "^0 

E / WlMfds < [ E ^m] ' [4(1 - X)-''Ge{Ue,V,) 
i=l "'O i=l 

+ 4X-\1 - X)-h^-Pgeiu„v,). (7.85) 

Proof. Dropping the indexes e,i for a moment and letting 0{s)'^ = 
(sina(s), — cosa{s))'^ we compute that 



and 



r \i{s)-^'{s)\ds < r XiM>l-X}{s)l'{s)-e{s)^ 

Jo Jo 

+2 [ X^M<l-x}{s) ds 
Jo 



: 2 - 2 det (y (5), 9{s)) = 2 sin (3{s) (—1^ - l) , 



ds 



(7.86) 



which imphes 



rX[M>i-x}{s) 7'{s)-e{s)^\s< 2(l-A)-2 ("^ (—}— i)M{sfds. 
Jo Jo ^sm/j^sj / 

(7.87) 



Moreover we calculate 

/ X{M<i-x}{s)ds < s^X-^ [ l(l-M(s))'ds. 
Jo Jo ^ 

By (7.87), (7.88) we obtain from (7.86) that 



(7.88) 



Jo 



|Y(s)-7'(s)|ds < \/2L£(l-A) 



1/2 



+e^X-^ I \{1-M{s)fds. 
Jo ^ 
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Summing the corresponding inequalities for i — 1, N{e), and using Holder's 
inequality, we deduce that 



Nis) 



i=l 



■{Jo £^ ^sin/3£,j(s) 



v)Me,i{sfds) 



1/2 



(7.89) 



(7.90) 



i=l 



where we have used (7.2) and (7.77). By (7.78) the inequality (7.83) follows from 
(7.90). 

Dropping indexes e, i again we obtain that 

["^irisWds = r XiM>i-x}{s)\a'{s)fds 
Jo Jo 

+ j X|M<,_j|(s)|a'(s)|>"<i« 

L Jo 4sm/3(s) \sm/3(s)/ 



+ 4(1- A)- 



eP 



'^{M<l-A}(s) ds, 



where we have used (7.80) in the last inequality. We can further estimate the second 
term on the right-hand side by (7.88), 

4(1 -A)-^^ t ^{M<i-x}{s)d8 < 4A-2(l-A)-2e2-f - M{s)f ds. 

s'^ Jo Jo ^ 

Putting everything together and summing over i = 1, N{e) we obtain 



E / ' ll'JM'ds < [E^ 



+ 4A-2(l-A)-V-f^,(w„z;,) 
which is (7.85). Setting p = 2 and using (7.2) the estimate (7.84) follows. □ 

The next result shows that those curves 7e,i that are 'too short' do not contribute 
to the limit. 

Lemma 7.21. There exists 6 = 6{A) and C = C(A) independent of e such that for 
the index set 



the estimate 



J{e) := {ie{l,...,7V(£)}:Le,i<<5} 

E Li^e < eC(A) 
ieJ{£) 



(7.91) 
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holds. 

Proof. Let £ > and i € J{s) be fixed and consider cig j and 'y^^i as in Definition 
7.18 with A = 1/2. We drop for a moment tlie indexes e, i. We tlien have 



L = 



fw{s)\^ds 
Jo 

I l\s)-{i{s)-m)' ds 
Jo 

- [ r{s) (lis) - 7(0)) ds + 7'(L) • (7(L) - 7(0)) 
Jo 



L 



L 



< / s\ns)\ds+ / (7'(,s) - y(s)) 
Jo Jo 

< ±L'/^(^£\^"{s)\Usy^\ £\j'{s)-j'{s)\ds. (7.92) 

Summing the corresponding inequalities over i G J{e) and using (7.83) and (7.84) 
we obtain with the Holder inequality 



£J(£) ieJ(e) 

< 6(^J2 ie,i)c(A) + eC(A) 
ieJis) 



ieJ(e) 

for S = S{A) small enough. This proves (7.91). □ 

Proposition 7.22. There exists a subsequence £ ^ and a W'^'^ -system T = 
{ci,i = 1, ...,Ni} of closed curves without transversal crossings such that 

s\\7ue\ Mr (7-93) 

as measures on . Moreover T has even multiplicity and satisfies 

W{T) < 2limm{g,{ue,v,). (7.94) 

e— >0 

Proof. By Lemma 7.21 we can, without loss of generality, assume that L^^i > S{A) 
for all e,i. We then deduce that the number of curves N{e) is bounded uniformly 
in £. Next we consider the reparametrized one-periodic functions Ce,, : M — » 
defined by 

Ce,j(r) 7e,i(L£,jr) for r G M. 

Since the functions c^^i are Lipschitz continuous with uniformly bounded Lipschitz 
constant, and since L^^i are bounded from above and away from zero, there exists 
a subsequence £ ^ 0, a number Ni gN and 



> 0, Cj : M ^ M^, i=l,...,Ni, 
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such that for all i G {1, A^i} the functions Ci are Lipschitz continuous and 



N{e) = Ni for all e > 0, 
Lg^i Li as e ^ 0, 

Ce,» Ci in C°''"(M;R2) as £ ^ 



(7.95) 
(7.96) 



holds for all < cT < 1. In particular, the functions Cj are one-periodic and elVuel 
converge as measures to /xr- 

We next turn to the VF^'^-bound on the limit curves We fix < A < 1 and 
consider for j^^i modified functions %^i as in Definition 7.18 and the reparametriza- 
tions Cs,i, 



leALe,ir) for r G R. 



By Lemma 7.20 those curves arc uniformly bounded in 14^j^'^^(M; R^) and we deduce 
that there exists a subsequence £ — > of the sequence in (7.96) such that the 
curves Ce^j converge weakly in W[^'^(M; M^) and strongly in C^''^ {[—K, KY.S?) for 
all < (7 < 1/2, K >Q. Since C£,i(0) = Cs,i{Q) we obtain for any r e M 

\im \ce,i{r) - Ce,i{r)\ < \im [ \c' i{g) - c' i{e)\ dg 



< hm f 1 + -P- 

£^0 V L^^i 

= 0, 



where we have used (7.83), the periodicity of Ce,i^ ^ and L^^i > (5(A). Therefore 
we can identify the hmits of c^^i, c^^i and deduce that 



Ci weakly in VKj^f ( 



as £ ^ 0, 



C^'''{[-K,K]-R^) as£^0. 



2,2/ 
loc ' 



for all < cr < 1/2, ii" > 0. In particular Cj € Wj 
semicontinuity of the norm under weak convergence and (7.85) yield that 

^Lj-^M c'lirfdr < hminf^^T^W c'^iirf dr 

i=l "'O ^ i=l "'O 



(7.97) 
(7.98) 

The lower- 



lim inf > 

£^0 ^ 



^1 rLe., 



^1 JO 



r ^1 nl-Er 

< '[4(1 -A)- 



lim \-niQs{ue,Ve) 

£— >0 



holds for all 1 < p < 2. Since < A < 1 and 1 < p < 2 are arbitrary we obtain 
(7.94). 

Since the curves Ce,i are pairwise disjoint and have no self-intersections, by (7.96) 
we deduce that the curves Ci have no transversal crossings. Finally we obtain for 
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any j] G C^(M^) by the divergence theorem that 
= hm / V • 'r]{x) dx 

J supp(w£ ) 



Ni .1 



Jsupp(r) L r,.- /_N 1 -'^i 



{{i,s):ci{s)=x} 



We deduce that 

E J^^'i^^^'^ = ^ for - almost all a; e supp(r) 

(i,s):c^(s)— a: * 

and since |c^| = Li and the vectors {c^(s) : Ci{s) = x} are collinear, this implies 
that 6r{x) = #{(«, s) : Ci(s) = x} is even. □ 

Using results from [13] we deduce that /ir as in Proposition 7.22 is given by an 
alternative system of curves, where each curve is passed twice. 

Lemma 7.23. Let T — {q, i = 1, .... TV} be a W'^''^ -system of closed curves without 
transversal crossings and with an even multiplicity function 6r- Then there exists 
a system of curves T = : i = 1,..., N} such that 

/ir = (7.99) 
W(T) = 2W(r), (7.100) 
6'r = 2%. (7.101) 

Proof. This Lemma requires some machinery from geometric measure theory. For 
the relevant definitions of Sobolev type submanifolds, Hutchinson varifolds and the 
various definitions for tangential lines and multiplicity functions see [13]. 
By [13] Remark 4.9 we obtain that 

/ := 6'(r) ■)'^supp(r) 

belongs to Sfg{M.'^), that is the set of Sobolev- type submanifold with the additional 
property that a unique tangential line exists in all points of the support. By [13, 
Proposition 4.5] this implies that / G HV^{M.^) is a Hutchinson varifold such that 
a unique tangential line exists in all points of the support. Since the multiplicity is 
even we conclude that also ^/ belongs to ffF^(M^) and that a unique tangential 
line exists in all points of the support. Using again [13, Proposition 4.5] we obtain 
that i/ e S^giR"^). Now [13, Corollary 4.8] gives that i/ is given by a VK^'^-system 
r = {7^ : i = 1, ...,N} of closed curves without transversal crossings. This implies 
(7.99). Equations (7.100) and (7.101) are an immediate consequence. □ 



We next show that the limit of the boundaries of the support of Ue and the limit 
of the mass distributions are identical. 
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Proposition 7.24. For a system of curves T as obtained in Proposition 7.22 and 
the measure n in (4.3) we have 

Mr = ^l■ 

In particular, (7.93) holds for the whole sequence e — > 0. 

Proof. Let 9supp(w,£) be given by a system of closed curves = {je,i}i=i,...,N(e) 
as in Remark 7.2. defined in Definition 7.3 we obtain from Lemma 7.10 that for all 

Z2 J \J v{le,i{^) + ts,i{^,m)ee,i{s))dm-ri{'ye,i{^))jds 



< 



\Me4s) - l\r]{-feAs)) ds 



E 



i=l JO JO 



??(7e,i(s) + te,i(s,TO)6'e,i(s)) - r]{'ye,i{^)) 

1 ,W(e) 





dmds 



<(X: / '"(M,,(s)-l)^d.)^(^ / "r,(7e,(.))'d.)^ 



+ ||'?||ci(R2) E / / te,i{s,m)dmds 
i=l -^0 ./o 

<£S£(W£,V£)2 2||?7||co(M2) + \\ll\\c^R^)dl{Ue,Ve), 

where we have used (7.71). We observe that by (7.2) 

and we deduce from (7.102) that 

/ vue- ??e|V«e| < £(2A^||7?||co(r2) + (2 + £2A)||77||ci(r2)) 



(7.102) 



This shows that 



/X = limuE^ = limelVuel = fir 

s—>0 e—^O 



in the sense of convergence as measures. 



□ 



The proof of the lim-inf estimate now follows: by Proposition 7.22 there exists 
a system of curves T of even multiplicity such that 

]-W{T)<limmig,{ue,Ve). 

Z £— >0 

By Lemma 7.23 there is an alternative system of curves F, of integer multiplicity, 

such that fir = 2/ip and yV'(r) = 2yV'(r). This second system F is the system 
referred to in Theorem 4.1 as F. Finally, Proposition 7.24 guarantees that /ir = 2/Ltp 
is the limit of the sequence {us,Vs) in the appropriate sense (see (4.3)). 
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8. The limsup estimate 

8.1. The heart of the construction. Let us first perforin the lim-sup construc- 
tion in the case that fi = 2/xr, where F consists of one simple curve. 

Lemma 8.1. Let 7 : [0, L) — !■ be a closed -curve without self-intersections 
and with \^'\ — 1 and let T — {7}. Then there exists eq > 0, £0 = £o(||7"||co([o,l))) 
and sequences (wg, We)o<e<eo 

e BF(M2;{0,l/e}), u,v,^Q, (8.1) 



/"' = / 



Ve = 2L, (8.2) 



supp(u£) U supp(w£) C {x : dist(x,r) < Se}, (8.3) 

such that 

2/ir = limue/;^ (8.4) 

£— >0 



Proof. There exists a Eq > such that 



\[ y (sf ds > limsup \(-di{u„Ve)+e [ \Vue\-iL). (8.5) 



sup \i\s)\ < (8.6) 

X dist(a;, T) is of class in {x : dist(a;, T) < 4eo}. (8.7) 
Let v : [0, L) — > be the unit normal field of 7 that satisfies 

det(7',i/) = -1 

and let k{s) be the curvature in direction of z^(s), 

k{s) = -u'{s)-j'{s) = iy{s)-j"{s). 

For £ < eo we set 

1 elsewise, 

and define two curves 7+, 7- : [0, L) ^ by 

7+(s) := 7(s)+£i/(s), 
7_(s) := 7(s)-£i^(s), 

see Figure 8. Using the parametrization $(s, r) := 7(5) + ri'{s) we calculate that 

I det Z)$(s, r) I = 1 - m(s) for < s < L, -£ < r < £. 

Therefore for any rj e C°(M2) 

lim / r}{x)ue{x)ds = lim / / ?7($(s,r))-(l — rK(s)) rfr ds 

= / 2ri{^{s,Q))ds = 2 1 r]{j{s))ds 
Jo Jo 

holds and we deduce (8.4). 
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The arc length functions corresponding to 7+ , 7- are given by 

= r\j'^{a)\da= f (l-£K(a))da, 
JO Jo 

PS PS 

^_(s) = / \f_{a)\da = / {l + SK{a))da. 
Jo Jo 

Let 7+,7_ be the corresponding parametrizations with respect to arclength, 

rL 



7+=7+o^;\ 7+:[0,L+)^M2^ L+ = L - e [ K{a) da, 

Jo 

7_=7_ot?l\ 7_ : [0,L_) IR^ L_ = L + s K{a)da. 

Jo 



The curves 7-|-,7_ parametrize the boundary of {u^ = l/s} and we obtain that 

e\Vue\ = L+ + L_ = 2L. (8.8) 



/ 



/M2 

Corresponding to 7+, 7- we define unit normal fields 

z>+ := i^oa?^^, 

:= v o {fZ^ 

and we compute the curvature of 7_|_ , 7_ in direction of 

K+ = — T> 

I — £K O 'tf_^ 

k- = t". 

1 + SKO 'dZ 

Observe that by (8.6) for e < eq 



\n+\,\K.\<^. (8.9) 
3eo 

As in subsection 7.1 we define mass coordinates m+,m_ by 

1 f 
m+{r,t) := -{t--k+{r)), 

1 

m_(r,i) := -(t--K_(r)) 

and the inverse mappings i+(r, •) = m+(r, •)~^, t-(r, •) = m_(r, ■)~^, 

t+(r,m) = . fl ~ ^1 - 2£k+m 
K,+ (r)\ 

t-(r,m) = — — 1 1 — ■\/T^^2sk~m) . 

K-{r) \ ) 

By (8.9) the expressions 1 — 2e'k±m are positive for |m| < 1. Using that 
Vl + 2z < 1 + ^, for2>-i, 



\/l + 2^ > 1 + 2^, for - - < ^ < 
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we deduce that 

< — t+(r,m) <2e, < — t_(r,m) < 2e. (8.10) 
m m 

We construct parametrizations ■;/'+, V"- iii the same way as we did in subsection 7.1. 
Considering 7+, the unit normal field !>+ takes the role of the direction field in 
subsection 7.1. Therefore k+ and 1 appear in place of a' and sin/3. Note that, 
with respect to 7+, supp('U£) is in the direction of — and, with respect to 7_, in 
direction of v- . 

For 7+,7'- respectively we define 'ray lengths' l+{r),l-{r) such that their values 




Figure 8. The limsup construction. 



in mass coordinates become —1 and 1, 

l+{r) :=t+(r,-l) = fl - ^1 + 2eR+) , (8.11) 

K+[r) \ / 

Z_(r) :=t_(r,l) = (l - ^1 - 2£k_) . (8.12) 

To parametrize supp(w£) we define sets 

L>+ := {(r, t) : < r < L+, /+(r) < t < 0} 
D- := {{r,t) : < r < L-,0 < t < l-{r)}. 

and the maps 

')p+{r,t) := 7+(r) +ti>+{r) for < r < L+,-2e < i < 2£: 

i)-{r,t) := 7-(r) +ti>_(r) for < r < L-,~2e <t <2e, 

which are, for e < Sq, by (8.7) injcctive on their domains. Analogous to Lemma 7.12 
we observe that for any 77 G '^0''»2^ 



i r]{x)us{x)dx = / r](^tp^{r,i^{r,m))) dmdr, (8.13) 

Ji)+(D+) Jo J-i 

I r]{x)ue{x)dx = / r]{^-{r,i-{r,m))) dmdr. (8-14) 

Jil>-{D-) Jo Jo 
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Moreover we deduce from (8.11), (8.12) that 



Z_ (^- («)) = ^ (l + ^«(«) - Vi-{eK{s)r) . 
which imphes that 

We conclude that the sets tp+{D+) and V'-(-D-) are pairwise disjoint and cover 
almost all of supp(u£). We therefore obtain by (8.13), (8.14) that 

ri{x)u^{x) dx 

= / / ri{^j^{r,i+{r,m))) dmdr + I / r][ip-{r,t-{r,m))) dmdr 
Jo J -I Jo Jo 

holds. In particular we obtain 

/ Ue = L+ + L_ = 2L. (8.15) 
Jm 

We now define an injective transport map S : supp(u£) — > IS? by 

S{tP+{r,t)) :=V+(r,t+(r,m+(r,i) + l)) for {r,t) G D+, 
5(V'-(r,t)) :=V-(r,t-(r,m-(r,t) - l)) for {r,t) e £>_ 

and set 

1 



^ '-S(supp(«e))- 



This gives 

supp{ve) = {ij+{r, t) : < t < i+{r, 1)} U {^'-(r, t) : t_(r, -1) < t < 0}. 

By (8.10) we obtain that t+(r, 1), |t_(r, — 1)| < 2e and since the distance of 7+ and 
7_ to 7 equals e the equation (8.3) follows. We compute for an arbitrary 77 e C^{M?) 
that 



y r]{S{x))ue{x)dx 



r](^S{ip+{r,t+{r,ni)))) dmdr + / / r][S{tp-{r,t-{r,m)))) dmdr 
-1 Jo -/o 



= I i ?7 (■!/)+ (r, t-|-(r,m + 1))) dm dr+ / / ?7(V'-(?', t_(r, m — 1))) rfmrfr 

Jo J-i Jo Jo 

= / r][ip-\.{r,t+{r,m))) dmdr + / / r][tp-{r,t-{r,m))) dmdr 

Jo Jo Jo J-1 

= J r]{x)vs{x)dx, 

where we have used the change-of-variables formula in the last equality. This 
shows that S indeed is a transport map. In particular it follows from (8.15) that 
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/r2 = 2L holds, which proves (8.2). Moreover we obtain 
di{us,Vs) < I \x - S{x)\us{x)dx 

= / (t+(r, m + 1) — t+(r, m)) dmdr 

+ / / {t-{r,m}-i-{r,m-l)))dmdr (8.16) 
Jo Jo 

By (8.11) we have m+(r, Z+(r)) = —1 and we deduce for the inner integral of the 
first term on the right-hand side of (8.16) that 



J (t+(r, m + 1) - t+(r, m)) dm 

/■m+(r,;+(r)) 

= / t+(r, m) — t+(r, m — m+(r, Z+(r))) rfm. 

Jo 



fm+(r,;+(r)) 
D 

Therefore Lemma 7.16 yields that 



.0 ^ 
/ (l+(r,m + 1) - t+(r, m)) dm = s + -e^k+{rf + R+ir)e^ 
J-i 4 



(8.17) 



where 

< R+{r) < K+(r)4 < C{eo) (8.18) 

by (7.73) and (8.6). For the inner integral of the second term in (8.16) we deduce 
by similar arguments 



1 

/ t_(r, m) — l_(r, to — 1) dm = e + -e'^k-{r)^ + R-{r)e^, 
Jo 4 



(8.19) 



with 

< R.-{r) < k^{rf < C(£o). (8.20) 
By (8.16) and (8.17)-(8.20) we obtain that 

^di{ue,Vs) 

<L+ + L_ + y^ ^ k+{r)^ dr + k-{rf dr + s'^LC{eo). (8.21) 
We compute that 



< 2 / K{sfds + 2e^LC{eo) 
Jo 



which gives with (8.21) and 2L = i+ + Z/_ that 

1 ,2 



-di{us,Vs) < 2L+ — K(s)^rfs + £UC(£o). (8.22) 
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Since \k\ = |7"| we obtain from (8.8) and (8.22) that 

^(-di{u„v,)+e [ |V«e|-4L) < i / ^"{sfds + e^LC{eo), 



which proves (8.5). □ 

8.2. The general case. We prove now the Hmsup estimate stated in Theorem 4.1. 
First we need a technical Lemma which follows from results in [11] and [12] and 
approximates a given system of curves as in Theorem 4.1 strongly in VF^'^ by a 
sequence of more regular systems of curves. 

Lemma 8.2. Let T be a W'^''^ -system of closed curves without transversal cross- 
ings. Then there exists a sequence (r-')jgN of W"^'^ -systems of closed curves and a 
number to S N such that the following holds: 

(a) For all j G N the system of curves — {7^}fc=i,...,m is given by a pairwise 
disjoint family of simple closed curves and satisfies 

|P'| = |r|. (8.23) 

(b) For all 1 < k < m holds 

Ik in W^'2(0, 1) as j ^ oo. (8.24) 

(c) T := {7fe}fc=i,...,m is a W^''^-system of closed curves without transversal cross- 
ings. 

(d) T is equivalent to T in the sense that 

Hr = Hf 

In particular we have 

W(P) ^ W(r) as j 00. (8.25) 

Proof. By [12, Corollary 5.2] there exists a sequence (r-')jgN and a W^^'^-system 
of closed curves F such that 

(1) for all j e N P is an oriented parametrization of a bounded open smooth set 

c IR^ 

(2) r-' converge weakly in W^'"^ to F as j ^ oo, 

(3) P converge 'in energy' to F as j — > oo. 

(4) F and F are equivalent. 

Property (1) implies that F' is a disjoint system of simple closed curves 
parametrized on (0, 1) proportional to arclength. 

By property (2) there exists a number to e N such that 

= {7fc}fe=l,...,m, f = {jk}k=l,...,m, 

7^,7fe:(0,l) ^ M2, 
\{%y\=Li |(7fc)'| =Lfeon(0,l), 

with 

7fe 7fe weakly in W'^''^{0, 1) as j oo. (8.26) 
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In particular we deduce that Lj, ^ Lfc as j ^ oo for all = 1, m and 

|P| - |f| (8.27) 

holds. 

The 'convergence in energy' stated in (3) gives that, as j ^ oo, 

|r^l+E / (LD-'Urr - |r| + E / (L0-'l(7fe)'f • (8.28) 

Together with (8.27) and the weak convergence (8.26) this yields 

7fe Ik strongly in W'^^'^{0, 1) as j oo, (8.29) 
>V(P) >V(r) as j oo. (8.30) 

Finally we define for j e N a modified system of curves = {7^}fe=i,...,m by 

Then (8.23) is satisfied. By (8.27) we see that 7^ and 7^ have for all 1 < fc < m the 
same limits as j 00. We therefore obtain from (8.29), (8.30) that (8.24), (8.25) 
hold. □ 



Proof of Theorem 4.1 (lim-sup part): 

Let F be given as in Theorem 4.1 and (r^)jgN be a sequence of systems of closed 
curves as in Lemma 8.2, given by F-^ = {7^}fe=i,...,m- We now parametrize the 
curves 7^ by arclength, 

7^(o,i:i)-M^ Hy\ = i. 

Choose Eg such that the distance functions dist(-,F-') is smooth in the set {x : 
dist(a;,rJ) < 44}- 

For < £ < £0 we let w^'*^, vi''' be the approximations constructed in Lemma 8.1 
and define 

m m 

fe=l fe=l 

By (8.3) all the functions have pairwise disjoint supports. We compute 

that 

/ uidC"^ ui^^dC^ = = 2|F^| = 2|F| = M 

"^■^^ fe=i "^^^ fe=i 

and, by similar calculations, j = M. Moreover, by (8.4), (8.24), 



m 
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and, by (8.25), (8.5), 
W(r) = lim W(P) 



1 _ _ r k 

> li^lsupVlimsup^(-rfl(^ij^^;p)+£ / |Vwi''=|-2 / vP^^ dC"^) 

> limsuplimsup ^ i;^) + £ / |Vw^| — 2 / uidC^\. 

j^OO E^O £ V£ J J 



3- 

The final inequality follows from constructing an admissible joint transport map 
for {u{,vi) from the transport maps of {ul''' , v^''') . 

Therefore there exist subsequences {ji)ieti , {£i)ien such that 

lim = 2/ir 

as Radon measures on and 

W(r) > limSUptJe, (Ue,,We,). 

I — >QO 

This proves the lim-sup estimate in Theorem 4.1. □ 

9. Discussion 

9.1. General remarks. In the Introduction we formulated two basic questions 
about self-aggregating, partially localized structures: first, why do they exist at all 
(and why are they stable), and secondly, how can we understand their resistance 
to stretching, bending, and fracture, that is observed experimentally. 

The fact that the functional J-'e favours structures in which u and v alternate on 
a length scale e is well illustrated by the special cases of Section 2. Simply put, for 
di{u,v) to be moderate, the masses represented by u and v have to be close; the 
fact that ||u||^oo and HfH^oo are limited implies that these masses occupy a certain 
amount of volume; and finally the interface penalization prevents the masses of 
u and V from being too finely interspersed. Together these restrictions impose a 
length scale on structures with moderate J^. 

This argument does not determine, however, the geometry of structures with 
moderate (or even lowest) energy. The fact that the straight lamellar geometry has 
lowest energy is not obvious, as is demonstrated by the 'wriggled lamellar' patterns 
of Ren and Wei in a strongly related system [42] . The reason for the penalization of 
curvature in the current context can be traced back to a property of strict convexity 
of the function t{s,m) defined in (7.54). Since this convexity property lies at the 
heart of the stability of these partially localized structures, we take some time to 
investigate it further. 

9.2. The geometric basis for stability. An essential observation in the proof 
is that it is more convenient to use mass coordinates along rays instead of length 
coordinates, since the Monge-Kantorovich distance measures spatial distances be- 
tween pairs of infinitesimal portions of mass. The function t(s, m) (see (7.54)) 
connects the two descriptions by giving the length coordinate t as a function of 
mass coordinate m. 
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For fixed m, the value of t(s,m) depends on the value of a'{s), as is obvious 
both from the definition and from a geometric point of view (see Fig. 9 (left)). 
For the purpose of this discussion we write the dependence on a'{s) explicitly, 
as t{s,m;a'{s)). The pertinent observation is that for fixed m 7^ 0, the function 
Q:'(.s) n- t(s, m; a'(s)) is strictly convex. This convexity can be seen in a simple plot 
of the function t, as in F'g. 9 (right), but can also be recognized in the geometry of 
rays that are rotated wil i respect to each other (Fig. 9 (left)). 



Figure 9. Left: the function t maps the mass coordinate m to 
a length coordinate t, in the following way. Fixing an infinitesi- 
mal section ds of interface, a section of mass mds, delimited by 
ds and by two rays, extends away from the interface by a dis- 
tance of t{s,m;a'{s)). Right: the function a'{s) 1-^ t(s,m;a'(s)) 
(see (7.54)), plotted for m G {1/3, 2/3, 1} (curves (1-3)). Here we 
have taken e = 1 and sin/3(s) = 1. 

With this convexity we can understand why curvature is penalized. In (7.75) 
the relevant expression is 



This integral describes the transport cost for the transport along an individual ray 
at position s. Using a symmetry property of t and the strict convexity mentioned 
above, 

t(s, m; Q;'(s)) — t(s, — m; Q;'(s)) = t(s, m; ct'{s)) + i{s, m; —a'{s)) 

> 2i{s,m; 0) 

= t(s, m; 0) — t(s, — m; 0). 

It is in this inequality that we see how curving the interface (or, more precisely, 
rotating the rays) creates a higher transport cost and therefore a larger value of the 
energy. 




o°.S 




which can be written as 




48 



MARK A. PELETIER AND MATTHIAS ROGER 



9.3. Connections to other work: continuum models and phase separa- 
tion. Beyond the direct context of lipid bilayers, this work is at the intersection of 
a variety of different hnes of research. There is of course a long-running tradition of 
slender-body limits in solid mechanics, resulting in a wide variety of simplified mod- 
els (sec e.g. [30] for an overview of the case of clastic plates and shells). Although 
our final result is rather similar — the Gamma-limiting energy W is the classical 
elastica bending energy — the current work is fundamentally different, in that the 
material at the e > level is closer to a fluid than a solid: particles may be ar- 
bitrarily rearranged without incurring an energy penalty. Indeed, the penalization 
of curvature that we find in the limit is an effect of global geometry rather than of 
(local) deformation. One might compare this with the difference between an elastic 
ball and a drop of water; both resist deformation, but for the ball this arises from 
local intermolecular forces, while the form of the drop results from the global effect 
of surface area minimization. 

Turning to fluid-like systems, the simple fact that the mass represented by u 
and V concentrates onto a curve is in itself remarkable. Phase separation naturally 
leads to penalization of interfacial length, as reflected by the term /|Vm| in J^g, 
and the distance penalization by the term di{u,v) is necessary to prevent bulk- 
scale separation. To our knowledge the current work is the first example of a 
phase-separating system that concentrates on low-dimensional sets. 

Concentration onto low-dimensional sets is a commonly observed feature in many 
(often non-conservative, non-variational) reaction-diffusion systems. For domains 
with sufficient symmetry a soliton- or pulse-like solution in one dimension can be 
interpreted as a solution in higher dimensions that concentrates on a hyperplane. 
In recent years less trivial examples of low-dimensional concentration have been 
uncovered [20, 9, 2, 3, 4, 23, 37, 38, 36]. Those concentrated solutions are often un- 
stable, however, and sometimes even highly unstable (e.g. [37, 38, 36]), thus strongly 
contrasting with the stable nature of the structures of this paper. 

9.4. Connections to other work: The elastica functional. Originally the 

elastica functional was introduced by Daniel Bernoulli and Euler as the bending 
energy of an elastic rod, but it has many other applications in different fields. In 
variational methods for image reconstruction it is widely used [39] . In the theory of 
phase separations the elastica functional appears as the two dimensional reduction 
of the Willmore functional (1.4). 

It is often natural to consider the elastica functional as acting on boundaries of 
subsets of the ambient space. In [11] the lower semicontinuity of the elastica func- 
tional under -convergence of sets is investigated; see also [44] for a generalization 
to arbitrary dimensions. The Gamma-convergence of a diffuse interface approxima- 
tion of the elastica, and more generally the Willmore functional, was conjectured 
by De Giorgi [22] and proved in a modified form in [43] for space dimensions two 
and three. 

Our approach is not restricted to curves which are boundaries of sets, and there- 
fore open curves can be represented in this framework. A related approach can be 
found in [15], where functionals which act on phase-field approximations of 'thick- 
ened curves' are considered and the Gamma-convergence of these functionals is 
proved in a topology based on the Haussdorf metric. The elastica functional is part 
of the limit, however the approximation is completely different. 
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9.5. Relevance to the understanding of lipid bilayers. The derivation that 
leads from a self-consistent mean-field theory for two-bead copolymers to the energy 
Te (Appendix A) contains a number of highly suspicious assumptions. The most 
glaring one is the severing of the head-tail bond in individual polymers: under this 
assumption a given head does not keep track of to which tail it is connected, as 
long as it is connected to exactly one tail. The Monge-Kantorovich distance d\(u, v) 
is the mathematical implementation of this assumption, and this term is all that 
keeps heads and tails from separating to arbitrarily large distances. 

The most interesting conclusion, from the point of view of applications, might 
be that this very weak remnant of the covalent bond suffices; that apparently very 
little is needed to convert a system exhibiting bulk-scale phase separation (e.g. an 
oil-water mixture) into a system which separates at a smaller scale. 

9.6. Other approaches to the upscaUng limit. In the physico-chemical liter- 
ature many authors have established connections between the Helfrich Hamilton- 
ian (1.3) on one hand and various smaller-scale models of lipid bilayers on the 
other [47, 41, 34]. Besides being an interesting scientific issue, such a connection 
should also give insight in the dependence of the coefficients k and k on molecular 
parameters; indeed the derivation given by Helfrich and others is phenomenological, 
based on scaling and dimensionality arguments, and therefore the coefficients are 
unknown. 

An approach that comes close to ours is that of [18], in which the authors connect 
the Helfrich description with a model in which lipids are represented by a density 

and a vectorial quantity that may be interpreted as an alignment order parameter. 
As in the examples above, however, the final coefficients are again found by a 
numerical method (minimization of a grand potential energy). Our result here 
appears to be; unique in giving a fully analytic formula for the bending modulus, 
which is not so surprising since the model contains no unscalable parameters other 
than £. 

Another way in which our contribution differs from existing work is the possibility 
to assess the energy penalties associated with non-optimal thickness and with curve 
'ends'. In the current paper we only show that these penalities are larger than 0(£^), 
but determining the actual scale of penalization is an interesting open question. 

9.7. Generalization to dp{u,v) for p > 1. The derivation of the energies 
presented in Appendix A suggests substituting the Monge-Kantorovich distance by 

the 2-Wasscrstein distance. Some of the results of this paper carry over to this 
situation, and in fact to the dp-metric for all 1 < p < oo. This follows from the 
remark that 



where we take total mass M equal to one for simplicity. Those statements of 
Theorem 4.1 that only depend on the upper bound (4.4) therefore carry over without 




= dp{u,v) 
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change for functionals such as 

J'P{u,v)=s j \\/u\ + ^dp{u,v) 

To be concrete, sequences {ue,Ve) along which the associated functionals or 
remain bounded will converge along subsequences to a system of closed VF^'^-curvcs. 

The Gamma-convergence result of this paper, however, can not be proved in this 
way: there is a gap between this lower bound and the upper bound that follows 
from a generalization of Section 8; this gap results from the strict convexity of the 
function z ^ ioT: p > 1. Nonetheless, we believe a similar result to be true for 
all 1 < p < DO. 

9.8. Generalizations. There are several other interesting directions in which fur- 
ther research could continue, among them 

• generalizations to higher dimension and higher codimension, 

• characterisation of local minimizers of the energies Te, 

• formulation of consistent evolution equations on the mesoscale, i.e. for the 
densities (u^, v^). 

In higher dimensions, and for limit structures of dimension larger than one, a con- 
cept similar to the 'systems of curves' will probably not be appropriate, and we 
expect that a varifold approach will be better suited. Some results of the present 
paper are easily generalized but there are also fundamental differences. 

To investigate how close the proposed mesoscale approximation is to the Helfrich 
energy or to the elastica functional it is interesting to compare (local) minimizers 
of the energy J^^ to the variety of shapes found as local minimizers for the Helfrich 
energy. 

A natural idea for extending the mesoscale description of lipid bilayers to a 
time-dependent evolution is to set up a gradient flow for the energies J-"e, such as 
in [28] or [14]. The latter corresponds to a gradient flow based on the 2-Wasserstein 
distance. 
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Appendix A. Derivation of 

The functional Te that is the basis of this paper arises in a simple model for a 
water-lipid system. In this model a lipid is represented by a two-bead chain: a head 
bead and a tail bead connected by a spring. The water molecules are represented 
by a third type of bead. 

The state space at the microscopic level is given by the positions X\, X"^, and 
of the lipid tail, lipid head, and water beads; the lipids are numbered by 
i = 1, . . . ,Ni and the water beads hy j = 1, . . . , Nyj. Assuming that the beads 
are confined to a space f2 c M'', the full microscopic state space for the system is 
then 

Elements X = {Xl , . . . , X^' , Xl, . . . , X^\ X^, . . . , X^-) G X are called mi- 
crostates. 

Wc describe the system in terms of probabilities on X, i.e. the (probabilistic) 
state %l) is & probability measure on X: 

ip&S, where f=|v:A'^M+, J ^ ^ 1^ . 

We assume that neither the microstates themselves nor the measure tjj can be 
observed at the continuum level; the observables are three derived quantities, the 
volume fractions of tails, of heads, and of water. 

For a given probability measure t/j G £, the water volume fraction rw{tlj) '■ ^ 
R+ is defined by 

ry,{tlj){x) = u V / t/; S{Xl - x) dX for all x G Q. 

Here v is the volume fraction of a single bead. The tail and head volume fractions 

r(('0) and rt{ip) arc defined similarly. 

To specify the behaviour of the system we introduce two free energy functionals. 
The 'ideal' free energy F"^ : £ ^ R models the effects of entropy and the interac- 
tions between beads of the same molecule; the 'non-ideal' free energy i^""^ : f ^ K 
represents the interactions between the beads of different molecules. The total free 
energy is the sum of the two, 

F{^) = F''^{^) + F'^''^{tp). 

For F"^ we assume zero temperature and postulate 

Jx 

where the function : A* ^ M is the internal energy of a microstatc. and the 
superscript 'id' again refers to a restriction to interaction within a single lipid 
molecule. While remarking that many different choices are found in the literature, 
here we choose simply to implement a spring by penalizing head-tail distance, 

k 

H'\x) = -Y,\xi-xi\r 

i=l 
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A natural choice is p = 2, which gives the energy of a hnear spring. Even more 
reahstic seems a general distance term that becomes infinite if \Xl — | exceeds 
a certain maximal value. 

For the non-ideal free energy F""^ we make an important simplifying assumption: 
F""^ can be written as a function of only the observables, 

Typical terms in the non-ideal energy are a convolution integral, in which 

proximity of hydrophilic (heads and water) beads and hydrophobic tail beads is 
penalized: 

/W / {r yj (ip) (x) +rh{ip)ix)) ri{ip){y)K{x - y)dxdy, 

and a compressibility term that penalizes deviation from unit total volume: 

^ / K(V')(a;) + n(V')(a;) + r^(V')(a;) - 1)^ dx. 
^ Jn 

In these expressions we take two limits: 

• In the limit K — > oo, we find that the observables satisfy an incompress- 

ibility condition 

+ rhW + r^ii;) = 1. (A.l) 
With this condition the convolution integral above can be written as 

M/ {l-rtW{x))rtiij){y)K{x-y)dxdy. (A.2) 
Jn Jn 

• Replacing the fixed function k by a rescaled version ks, defined by 

Ks(x) :— S~'^k(x/S), 
the integral (A.2) Gamma-converges in the following way, 

n(V')(y) i^s{x - y) dxdy ^ (^j F'"*(n(V)) 

where the limit interface functional F™^ is given by 

1 oo otherwise. 

Here Ki^i = 1, Ki^2 = 2/7r, and i^i,3 = 1/2 [1, 21]. 
Putting the ideal and non-ideal energies together, and applying the simplifica- 
tions above, we find (up to rescaling) 

F{i,) := + F^\i,) = / ^|x;-X^|^V(^)dX + F-'(n(V')). 

We now assume that we minimize F under all variations of tj) that conserve the 

observables riji_u,{ijj). By a convexity argument, using the fact that the fimctions 
f't,h,w do not distinguish between different lipids or water molecules, it follows that 
all lipids are identically distributed, as are all water molecules: 

^{X) = l[MXl,Xi,)l[i;^{Xi,), 
i=i j=i 
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and that the total energy becomes a function of the observablcs alone, by mini- 
mization over all other degrees of freedom (again we disregard rescaling) 

F{puPh,Pw) ■■= F"'\pt) + 

+ inf I / \Xt - Xhf MXt, Xh) dXtdXh : = pt, rh{i^e) = Ph] 

with the remaining constraint pt + ph + Pw = ^- Remark that the final infimum is 
exactly the Wasserstein distance of degree p, dp{pt,phY- 

We now recover the functional Te in (1.1) for e = 1 by defining u := pt, v := ph 
and by choosing p = 1. Combining the pure-phase condition u{x) — pt{x) G {0, 1} 
for all X (A. 3) with the incompressbility condition (A.l) gives v{x) = Ph{x) € [0, 1] 
and uv = Q. It is easy to see that there is nothing to be gained in letting v take 
values in the interior (0, 1), and we can therefore consider to be given on the set 
of admissible functions 

Ki := ^{u,v) e BV(K^{0,1}) x L^{^^;{0,1}) : 

j ^ — J V = M, uv = a.e. I . 
The remaining dependence on e is a consequence of simple rescaling. 

Appendix B. Ring solutions 
We consider ri < r2 < R < < r4 with 

^{R'+rf), rl = UR^ + rl), (B.l) 



r. 



2R = ri+ r4 (B.2) 

and u, V given by 

In this setting the unique monotone optimal transport map S from u to t; is radially 
symmetric, S{x) = S{\x\) and determined by the requirement that 

) 

2'Krdr for < s < i?, (B.3) 
2irrdr for ii < s < ra. (B.4) 

Js JS{s) 

Then 

di{u,v) = J \x — S{x)\u{x) dx 

j-R j-rs 

= / 2'Kr{r - S{r))dr + I 2'Kr{S{r) - r) dr. (B.5) 

Jr2 Jr 

To compute the first integral we introduce new coordinates 

m(r) = -KT^ —irr^ with inverse r{m) = (r2 H — m)^^^. (B-6) 





rS(s, 


/ 2nrdr = 




J r2 


Jri 






' 2'Krdr 




■ 


Js{s) 
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In this coordinates (B.3) implies 

m(r) = m{S{r)) — m(ri), 
S{r{m)) = r[m + m{ri)) = r[m — m{R)) , 

where in the last equality we have used that 

m(ri) + m{R) = 7r(r? - 2r^ + R^) =0 

by (B.l). The first integral in (B.5) therefore becomes 

/ 2irr{r — S{r)) dr = (r(m) — r{m — m{R))) dm 

Jr2 Jo 



(B.7) 
(B.8) 



I 



m{R) 



r{m) — r{—m) dm 





27r r 

Y 
Y 



2^ 1 

ro H — m 



3/2 



3/2- 



V 1 / 1 \3/2 

i^rl + -m{R)j + i^rl - -m{R)j - 2r: 
If we now introduce t = (r4 — ri)/2 then 



(B.9) 



R — Ti = r4 — R = t, 



m{R) = '■^{2Rt-t^), 



and (B.9) yields 



/' 



27r 



27rr(r - S{r)) dr =^ R^ + {R- tf - 2R: 



y-R^2^) 



(B.IO) 



If we do the analogous calculations for the second integral in (B.5) we obtain that 



JR 



27r r 

2TTr{S{r) -r)dr = — R^ + {R + tf - 2R^ 



{'+R + 2¥) 



(B.ll) 



which is the expression on the right-hand side of (B.IO) with t substituted by —t. 
We consider R oo with t of order one and obtain by a Taylor expansion in 1 /R 

di{u,v) + j |Vw| 



27r 

Y 



+ 2-KR 



t_ 



i?^2i?2; 



= 2tt Rt 

= 2TTRt 



t 2 



4^ + W + 2* - 4) + 0(i?-^)0(|t - 2|) 



2 i 

2 + i(i - 2f + ]r-^ + 0{R-^)0{\t - 2|) 



4' ' 4 
which is the desired asymptotic development. 



partial localization, lipid bilayers, and the elastica functional 55 

References 

[1] G. Alberti, G. Bellettini, M. Cassandro, and E Presutti. Surface tension in Ising systems with 
Kac potentials. J. Stat. Phys., 82:743-796, 1996. 

[2] A. Ambrosetti, A. Malchiodi, and W.-M. Ni. Solutions, concentrating on spheres, to sym- 
metric singularly perturbed problems. C. R. Acad. Sci. Paris, Ser. I, 335:145-150, 2002. 

[3] A. Ambrosetti, A. Malchiodi, and W.-M. Ni. Singularly perturbed elliptic equations with 
symmetry: Existence of solutions concentrating on spheres, Part I. Commun. Math. Phys., 
235:427-466, 2003. 

[4] A. Ambrosetti, A. Malchiodi, and W.-M. Ni. Singularly perturbed elliptic equations with 
symmetry: Existence of solutions concentrating on spheres. Part II. Indiana Univ. Math. J., 
53:297-329, 2004. 

[5] L. Ambrosio. Lecture notes on optimal transport problems. In Mathematical aspects of evolv- 
ing interfaces (Funchal, 2000), volume 1812 of Lecture Notes in Math., pages 1-52. Springer, 
Berlin, 2003. 

[6] L. Ambrosio. Optimal transportation and applications, volume 1813 of Lecture Notes in Math- 
ematics. Springer- Verlag, Berlin, 2003. Lectures from the C.I.M.E. Summer School held in 
Martina Franca, September 2—8, 2001, Edited by Caffarelli and S. Salsa. 

[7] L. Ambrosio, N. Fusco, and D. Pallara. Functions of Bounded Variaton and Free Disconti- 
nuity Problems. Oxford Science Publications, 2000. 

[8] L. Ambrosio, N. Gigli, and G. Savare. Gradient flows in metric spaces and in the space of 
probability measures. Lectures in Mathematics ETH Zurich. Birkhauser Verlag, Basel, 2005. 

[9] M. Badiale and T. D'Aprile. Concentration around a sphere for a singularly perturbed 
Schrodingcr equation. Nonl. Anal, 49:947-985, 2002. 
[10] F. S. Bates and G. H. Frcdrickson. Block copolymer thermodynamics: Theory and experi- 
ment. Annu. Rev. Phys. Chern., 41:525-557, 1990. 
[11] G. Bellettini, G. Dal Maso, and M. Paolini. Semicontinuity and relaxation properties of a 
curvature depending functional in 2D. Ann. Scuola Norm. Sup. Pisa CI. Sci. (4), 20(2):247- 
297, 1993. 

[12] G. Bellettini and L. Mugnai. Characterization and representation of the lower semicontinuous 
envelope of the elastica functional. Ann. Inst. H. Poincare Anal. Non Lineaire, 21(6):839- 
880, 2004. 

[13] G. Bellettini and L. Mugnai. On the approximation of the elastica functional in radial sym- 
metry. Calc. Var. Partial Differential Equations, 24(1):1— 20, 2005. 
[14] J. G. Blom and M. A. Peletier. A continuum model of lipid bilayers. Euro. Jnl. Appl. Math., 

15:487-508, 2004. 

[15] A. Braidcs and R. March. Approximation by F-convergence of a curvature-depending func- 
tional in visual reconstruction. Comm. Pure Appl. Math., 59(1):71-121, 2006. 

[16] Luis A. Caffarelli, Mikhail Feldman, and Robert J. McCann. Constructing optimal maps for 
Monge's transport problem as a limit of strictly convex costs. J. Amer. Math. Soc, 15(l):l-26 
(electronic), 2002. 

[17] P.B. Canham. The minimum energy of bending as a possible explanation of the biconcave 

shape of the human red blood cell. J. Theor. Biol., page 61, 1970. 
[18] E. Chacon, A.M. Somoza, and P. Tarazona. Elastic constants from a microscopic model of 

bilayer membrane. J. Chem. Phys., 109:2371-2379, 1998. 
[19] B. Dacorogna and W. Gangbo. Extension theorems for vector valued maps. 

www.math.gatech.edu/~gangbo/publications/extension.pdf. 
[20] T. D'Aprile. Behaviour of symmetric solutions of a nonlinear elliptic field equation in the 

semi-classical limit: Concentration around a circle. Elec. J. Diff. Eqns., 69:1-40, 2000. 
[21] J. Davila. On an open question about functions of bounded variation. Calc. Var., 2002. 

published electronically. 

[22] E. De Giorgi. Some remarks on F-convcrgcncc and least squares method. In Composite me- 
dia and homogenization theory (Trieste, 1990), volume 5 of Progr. Nonlinear Differential 
Equations Appl., pages 135—142. Birkhauser Boston, 1991. 

[23] A. Doelman and H. van der Ploeg. Homoclinic stripe patterns. SIAM J. Appl. Dyn. Sys., 
1:65-104, 2002. 

[24] E. Evans. Bending resistance and chemically induced moments in membrane bilayers. Bio- 
phys. J., 14:923, 1974. 



56 MARK A. PELETIER AND MATTHIAS ROGER 

[25] L. C. Evans and W. Gangbo. Differential equations methods for the Monge-Kantorovich mass 
transfer problem. Mem. Amer. Math. Soc, 137(653) :viii+66, 1999. 

[26] L. C. Evans and R. F. Garicpy. Measure Theory and Fine Properties of Functions. CRC 
Press, Boca Raton - Ann Arbor - London, 1992. 

[27] M. Feldman and R. J. McCann. Uniqueness and transport density in Monge's mass trans- 
portation problem. Calc. Var. Partial Differential Equations, 15(1):81-113, 2002. 

[28] J. G. E. M. Praaije. Dynamic density functional theory for microphase separation kinetics of 
block copolymer melts. J. Chem. Phys., 99:9202-9212, 1993. 

[29] G. H. Predrickson and F. S. Bates. Dynamics of block copolymers: Theory and experiment. 
Annu. Rev. Mater. Sci., 26:501-550, 1996. 

[30] G. Friesccke, R. D. James, and S. Miillcr. A hierarchy of plate models derived from non- 
linear elasticity by gamma-convergence. Technical Report 7/2005, Max Planck Institute for 
Mathematics in the Sciences, Leipzig, 2005. 

[31] W. Helfrich. Elastic properties of lipid bilayers: Theory and possible experiments. Z. Natur- 
forsch. Teil C, 28:693-703, 1973. 

[32] R. Jordan, D. Kinderlehrer, and F. Otto. The variational formulation of the Fokker-Planck 
equation. SIAM J. Math. Anal, 29(1):1-17 (electronic), 1998. 

[33] L. V. Kantorovich. On a problem of Monge. Zap. Nauchn. Sem. S.-Peterburg. Otdel. Mat. 
Inst. Steklov. (POMI), 312(Teor. Predst. Din. Sist. Komb. i Algoritm. Metody. 11):15-16, 
2004. Reprinted from C. R. (Doklady) Acad. Sci. URSS (N.S.) 3 (1948), no. 2. 

[34] M. Laradji and O.G. Mouritsen. Elastic properties of surfactant monolayers at liquid-liquid 
interfaces: A molecular dynamics study. J. Chem. Phys., 112:8621-8630, 2000. 

[35] L. Leibler. Theory of microphase separation in block copolymers. Macromolecules, 6:1602— 
1617, 1980. 

[36] A. Malchiodi. Concentration at curves for a singularly perturbed neumann problem in three- 
dimensional dom ains. Geom. Funct. Anal., 15(6):1162-1222, 2005. 

[37] A. Malchiodi and M. Montenegro. Boundary concentration phenomena for a singularly per- 
turbed elliptic problem. Comm. Pure Appl. Math., 55:1507-1568, 2002. 

[38] A. Malchiodi and M. Montenegro. Multidimensional boundary layers for a singularly per- 
turbed Neumann problem. Duke Math. .J., 124(1):105-143, 2004. 

[39] D. Mumford. Elastica and computer vision. In Algebraic geometry and its applications (West 
Lafayette, IN, 1990), pages 491-506. Springer, New York, 1994. 

[40] F. Otto. The geometry of dissipative evolution equations: the porous medium equation. 
Comm. Partial Differential Equations, 26(1-2): 101-174, 2001. 

[41] S.M. Oversteegen and E.M. Blokhuis. Rigidity constants from mean-field models. J. Chem. 
Phys., 112:2980-2986, 2000. 

[42] X. Ren and J. Wei. Wriggled lamellar solutions and their stability in the diblock copolymer 
problem. SIAM I. Math. Anal., 37(2):455-489 (electronic), 2005. 

[43] M. Roger and R. Schatzle. On a modified conjecture of de giorgi, 2006. submitted. 

[44] R. Schatzle. Lower semicontinuity of the willmore functional for currents, 2004. submitted. 

[45] U. Seifert. Configurations of fluid membranes and vesicles. Advances in Physics, 46(13), 1997. 

[46] L. Simon. Lectures on Geometric Measure Theory. Proceedings of the Centre for Mathemat- 
ical Analysis Australian National University Vol. 3, 1983. 

[47] I. Szleifer, D. Kramer, A. Ben-Shaul, W.M. Gelbart, and S.A.Safran. Molecular theory of 
curvature elasticity in surfactant films. J. Chem. Phys., 92:6800-6817, 1990. 

[48] N. S. Trudinger and X.-J. Wang. On the Monge mass transfer problem. Calc. Var. Partial 
Differential Equations, 13(1):19-31, 2001. 

[49] C. Truesdell. The influence of elasticity on analysis: the classic heritage. Bull. Amer. Math. 
Soc. (N.S.), 9(3):293-310, 1983. 

[50] C. Villani. Topics in optimal transportation, volume 58 of Graduate Studies in Mathematics. 
American Mathematical Society, Providence, RI, 2003. 

[51] T.J. Willmore. Note on embedded surfaces. Ann. Stiint. Univ. Al. I. Cuza, lagi, Sect. I a 
Mat. (N.S.), llB:493-496, 1965. 



